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ABSTRACT 


The flowfield about a vertically-launched surface-to-air 
missile model at an angle of attack of 50 degrees and a 
Reynolds number of 1.1x10° was investigated in a low-speed 
Wind tunnel at the Naval Postgraduate School. The location 
and intensity of the asymmetric vortices in the wake of the 
missile model were determined and the vortices were displayed 
using planar velocity vector, total pressure coefficient, and 
vorticity plots. The model configuration tested was a body- 
only configuration (wings, strakes, and tails removed). Two 
flowfield conditions were treated: the nominal ambient wind 
tunnel condition and a grid-generated turbulence condition. 
Flow visualization was conducted and video-taped for both the 
body-only configuration and the winged configurations. The 
following conclusions were reached: 1) the addition of 
turbulence decreased the vorticity but did not significantly 
change the patterns of the plots; 2) the addition of 
turbulence reduced the vorticity more at eleven body diameters 
than at six body diameters; 3) compared to the body-only case, 
the vorticity is reduced for the "x" case but not for the "+" 
case for the turbulence condition; 4) flow visualization 
verified vortices movement away from the missile as the tested 


point was moved aft along the missile body. 
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I. INTRODUCTION 


A. BACKGROUND 

Missile development in recent years has emphasized the 
capability for operation at high angle of attack due to demand 
for higher maneuverability and special mission requirements. 
This trend is amplified by the increasing demand to minimize 
the overall wingspan. Several examples of missile design 
which involve high angles of attack are vertical launching, 
missile storage separation from aircraft, launch from an 
aircraft container, the one-man anti-tank weapon, and the air- 
to-air and air-to-ground short range missile. These are 
examples where high angles of attack can be reached at 
subsonic Mach numbers. The flight Mach number in some cases 
can even range from low subsonic to high supersonic. In such 
flight conditions it is important to understand the flowfield 
around the missile and the effects on aerodynamic 
characteristics. [Ref. 1] 

An example of such a missile is the U.S. Navy's 
vertically-launched surface-to-air missile (VLSAM). Upon 
launch it enters the ocean's surface environment at low 
velocity where it may encounter significant crosswinds 
resulting ina high angle of attack situation. In addition, 


the VLSAM launch conditions may include some degree of 


turbulence caused by interaction of the air movement above the 
ocean surface, the ship's = Metadonm and superstructure 
interference. [Ref. 2] High angle-of-attack conditions may 
result in the formation of asymmetric vortices along the 
missile body. These vortices induce side forces on the 
missile and can affect its flight aerodynamics. A vertically- 
launched anti-submarine rocket was lost due to asymmetric 
vortices caused by its high angle of attack and nose asymmetry 
ker. 34). 

This thesis is a part of the research conducted at the 
Naval Postgraduate School (NPS) to understand the effects of 
turbulence on a VLSAM model at high angle of attack. 
Initially, Roane developed the missile model and procedure to 
measure various flowfield turbulence effects [Ref. 4}. Rabang 
studied the asymmetric-vortex-induced side force by varying 
turbulence, angle of attack, and nolieanai. He concluded 
that turbulence with length scales comparable to the nose- 
generated vortices increased the side force and decreased flow 
unsteadiness. He also verified that nose roll angle can cause 
Significant changes in the side forces and that the addition 
of low-aspect wings did not significantly change the magnitude 
of the side forces. [Ref. 2] 

Lung mapped the location and shape of the asymmetric 
vortices, in a cross-plane six missile diameters from the 


nose, using a body-only configuration with and without 


freestream turbulence, at 50 degree angle of attack. He 
concluded that the strength of the vortices weakened and 
became more diffuse with turbulence. [Ref. 5] 

Viniotis repeated the procedure at six diameters using two 
wing configurations: zero degree Old angle (uy 
configuration) and 45 degree roll angle ("x" configuration). 

He concluded that the addition of turbulence decreased vortex 
strength and the addition of wings causes the vortices to move 
closer to the missile body. [Ref. 6] 

Johnson continued the study by mapping the flowfield in 
a cross-plane eleven diameters from the nose for the "+" and 
oe" configurations. He concluded that the addition of 
turbulence decreased the vortex strength with little change 
in the vortices' positions. He also noted that the vortex 
strength was greater for the "+" configuration than for the 
"x" configuration, but that the amount of change was much less 
with turbulence addition. [Ref. 7] 

The purpose of this research was to study the vortices for 
the body-only configuration in turbulent and non-turbulent 
flow at eleven missile diameters from the missile nose. [In 
addition, flow visualization was conducted for all missile 
configurations. Comparisons with the results of Lung and 
Johnson will be made io) correlate the vortices' 
characteristics at six diameters and eleven diameters for all 


missile configurations. 


B. HIGH ANGLE OF ATTACK AERODYNAMICS 

An important phenomenon of high angle-of-attack 
aerodynamics is the flow separation on body, wing, and tail 
surfaces and the shedding of free vortices. The flow 
Separates along a line on the body surface via an interaction 
of the external flow with the laminar or turbulent boundary 
layer of the body. The separation causes the formation of 
feeding sheets which transfer the vorticity generated in the 
boundary layer into the external flow, which then roll up on 
the lee side of the body to form vortices. [Refs. 1,8] 

A slender body of revolution experiences four distinct 
aerodynamic regimes as the angle of attack is increased from 
O to 90 degrees. The transition from one regime to the next 
as a function of angle of attack is also dependent on nose 
Shape, overall fineness ratio, crossflow Mach number, Reynolds 
number, roll angle, free stream turbulence, and surface 
roughness. [{Refs. 9,10] 

Figure 1 shows schematics of the four flow patterns about 
an ogive-cylinder body and the corresponding characteristics 
of normal force and side force. At very low angles of attack 
(0 to 5 degrees) there is no noticeable boundary layer 
Separation. In this regime the axial flow components dominate 
and the flow can be characterized as a classical potential 
flow field and an attached laminar or turbulent boundary 


layer. At intermediate angles of attack (5 to 20 degrees) 
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Figure 1. Flow Regimes ([Ref. 2] 


the boundary layer separates on the lee side of the body and 
rolls up into a symmetric vortex pair that originates from the 
apex of the body and extends to the base of the body. The 
vortices are steady with time and no side force or yawing 
moment is present. [Ref. 9] 

In addition to the primary symmetric vortex pair, a second 
pair of vortices exists with the same rotational direction as 
the primary pair, positioned near the separation line as shown 
in Figure 2. The second interior vortices pair and the 
secondary vortices are induced by the primary vortices with 
rotation directions opposite to those of the primary vortices. 
The primary and the second vortices move away from the body 
and their cores become larger as the angle of attack 
increases. The second and secondary vortices are usually 
negligible in strength as compared to the primary vortices. 
PRetedia | 

In the next regime, 20 to 60 degrees, crossflow effects 
begin to dominate and the vortices may become asymmetric, 
producing a side force and yawing moment. The formation of 
asymmetric vortices is discussed in the next section. Lastly, 
at an angle of attack between 60 and 90 degrees, the crossflow 
completely dominates and the vortex shedding becomes 
unsteady, starting aft and moving forward as the angle of 
attack increases. The boundary layer is shed in the form of 


a Karman vortex street or random wake depending upon the 


second vortex 


= 
Primary VOrtox l{A/ 
( we, 4] secondary 


vortex 






Figure 2. Cross~-flow Streamlines [{Ref. 11] 


Reynolds number, Mach number, body geometry, etc. [{Refs. 


SP 10,11) 


C. ASYMMETRIC VORTICES 

As the angle of attack transitions through the 20 to 60 
degree regime, the vortex flowfield on the lee side of slender 
bodies may become asymmetric, even without sideslip. These 
asymmetric vortices produce an undesirable side force and 
yawing moment which can make necessary a complex autopilot or 
degrade the mission effectiveness of the missile. The side 


forces may cause an unacceptable target miss distance, and the 


yawing moments may reduce stability and controllability. 
[(Refs.9,11] Reynolds number and Mach number affect the 
magnitude of the side force, but do not strongly influence the 
angle of attack for the onset of the side force. In general 
an increase in Mach number will decrease the magnitude of the 
Side force. [Ref. 12,13] 

The cause of asymmetric vortices formation is. not 
completely understood, but there are two generally accepted 
theories. The first theory is that boundary-layer-induced 
asymmetry at the flow separation points causes the vortex 
flowfield to become asymmetric. The maximum vortex-induced 
normalized side force will occur in the critical Reynolds 
number regime where it is possible to have supercritial 
separation on one side of the body and subcritical separation 
on the other. [Refs. 14,6] The second theory suggests the 
principle cause of vortex flow asymmetry to be a hydrodynamic 
(inviscid) instability in the initially symmetric vortex 
formation. These vortices, which increase in strength with 
angle of attack, interact with the surrounding potential 
flowfield to form the asymmetric characteristic. [Refs. 12,15] 

A vortex-switching phenomenon has been observed in which 
the vortex pattern rapidly switches from an almost symmetric 
to a highly asymmetric configuration [Ref. 13]. This 
Switching phenomenon may result from a fluctuation in the 


separation characteristics on the missile, or may relate to 


a second inviscid solution in the leeward flowfield [Ref. 16]. 

The large "steady" side forces observed by many investigators 
indicate that the detached vortices fluctuate about some mean 
asymmetric configuration [Ref. 13]. 

The missile fore-body has proven to be the dominating 
component concerning flow separation [Ref. 1]. The addition 
of wings does, however, cause the vortices to move closer to 
the body. This results in an increase of the effective angle 
of attack for the formation of the asymmetric vortices. [Ref. 
18] The addition of strakes to a missile body greatly 
increases the lift and generates additional vortices which may 
alter the effect of the asymmetric vortices. The addition of 
tails has little effect on the flowfield and side forces of 
the forebody. Rabang has shown that the missile model retains 
the side force and length scale effects (described in the next 
section) of the induced vortices with the addition of strakes 


and wings. [Ref. 2] 


D. TURBULENCE 

Turbulence represents the presence of random, short 
duration variations ina flowfield with a given mean velocity. 
When considering the effects of turbulence on a body in a 
flowfield, it is important to compare the scale of the body 


to that of the turbulence. [Ref. 2] 


Turbulence intensity measures the relative magnitude of 
high frequency velocity fluctuations in the flowfield. 
Greater intensities correspond to more kinetic energy in the 
turbulent § flow. The turbulence intensity, 2) 
mathematically expressed by the ratio of the streamwise root- 
mean-square (rms) velocity fluctuation, u', to the mean 
velocity component of the flowfield, Um. [Ref. 2] 

Ti =u! /Ue (1) 

Turbulence length scales are a measure of the time- 
averaged size of the fluid disturbance eddies. The turbulence 
length scale to body size ratio may determine the amount to 
which the turbulence affects the missile in the flowfield. 
Comparing the length scale to the missile length or missile 
diameter gives three possible results. For a length scale 
much greater than the body, the effect on the missile is 
Similar to that of a steady-state flowfield. For a length 
scale on the same order as the body, the flowfield is non- 
steady and may cause rolling, pitching and yawing motion of 
the body. [Ref. 2] For a length scale much smaller than the 
body, most noticeably smaller than the missile diameter, the 
small-scale turbulence may affect the boundary layer 
development and flow separation over the missile body. [Ref. 
19] 

Since the length scale is an average of the disturbances, 


turbulence length scales of varying magnitude can be found in 


me 


a flowfield. These magnitude differences can be explained by 
the cascade effect. The turbulence eddies in the flowfield 
experience a strain in one direction which affects the scale 
of motion in the other directions as required by the 
conservation of angular momentun. The strain causes the 
eddies to break up into smaller scale and decreasing energy 
disturbances. The "cascade" of energy of turbulent motion 
continues to smaller and smaller scales until viscosity 
dissipates the smallest eddies. As the turbulence decreases 
in size and energy, the individual intensities decrease at a 
faster rate. These changes account for the dominance of the 


larger scale turbulence. [Ref. 20] 
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II. EXPERIMENT AND PROCEDURES 


A. OVERVIEW 

Pressure measurements were recorded in the VLSAM model 
wake to study the behavior of the asymmetric vortices. The 
body-only missile model was positioned at 50 degrees angle of 
attack for freestream conditions with and without generated 
turbulence in the Naval Postgraduate School low-speed wind 
tunnel. The pressure measurements were converted to total and 
planar velocity components, total pressure coefficients, and 
vorticity which were plotted to observe the vortices' 
characteristics. 

Figure 3 shows the survey grid in an x-y plane which was 
perpendicular to the freestream flow 11 missile diameters aft 
from the missile nose, measured along the missile body. The 
Survey grid measured 3.5 inches by 5.0 inches. Incremental 
step size was 0.25 inches, yielding 15 columns in the 
horizontal direction each having 21 points in the vertical 
direction. The grid was centered on the missile body and 
started 1.50 inches away from the missile body due to physical 


constraints of the sting mounting arm and probe. 
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(b) Front View 





Planar Survey Grid [Ref. 5] 
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B. APPARATUS 
Primary equipment used were the NPS wind tunnel, 
turbulence-generating grid, VLSAM model, three-axis traverser, 
five-hole pressure probe, data acquisition system, and data 
reduction/display software. 
1. Wind Tunnel and Turbulence Grids 
The low-speed, single-return, horizontal-flow wind 


tunnel is located in Halligan Hall at NPS (Figure 4). 
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Figure 4. Naval Postgraduate School Wind Tunnel [Ref. 21] 
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The tunnel is powered by a 100-horsepower electric 
motor with a three-blade variable pitch fan and a four-speed 
Dodge truck transmission. Stator blades aft of the fan help 
to straighten the flow. Turbulence is reduced by two wire 
mesh screens upstream of the settling chamber and corner 
turning vanes. The test section is 45 inches by 32 inches 
with corner lighting enclosures that reduce the section area 
from 10 ft* to 9.88 ft*. A reflection plane in the test 
section decreases the useable height to 28 inches. The 
tunnel's contraction ratio is approximately 10:1. The tunnel 
was designed for test section velocities up to 290 ft/sec, 
operating at atmospheric pressure. A remote-controlled, 
flush-mounted turntable in the test section allows changing 
of model angle of attack in a horizontal plane. Wind tunnel 
temperature was measured in the settling chamber with a dial 
thermometer. Dynamic pressure was determined by the pressure 
difference between the static pressure in the settling chamber 
and in the test section by a water micro-manometer. The 
settling chamber and test section each have four static taps 
connecting to the manometer via a common manifold. [Ref.21]} 
The measured pressure difference, measured in centimeters of 
water by the manometer, is converted to the test-section 
reference velocity by equation (2). 


Uy=( (2*2.046*cm H,0)/(K*p) ) |” (a 


1: 


Where: 
Oe = measured velocity (ft/sec) 
2.046 = conversion factor 


cm H,O = manometer reading in cm of H,O 


K = wind tunnel calibration factor (for specific 
epaayel) 
0 = air density (loses) 


The wind tunnel calibration factor, K, corrects@iias 
the fact that the actual dynamic pressure is slightly 
different than the measured static pressure difference between 
the static pressure in the settling chamber and in the test 
section. The value of the calibration factor was determined 
by plotting the actual dynamic pressure measured by a pitot- 
static tube in the test section versus the measured pressure 
difference. (Ref. 21] The relationship was found to be linear 
where the slope of the curve is the tunnel calibration factor. 
The value of K was determined to be 0.8891 for the no-grid 
configuration and 1.6545 for grid three. A large wall 
disturbance due to the grid frame accounts for the change in 
Ka) (Ret. Gi 

Four square-mesh, square-bar, biplanar grids were 
designed and built by Roane to generate turbulence of varying 
intensities and length scales. Each grid, three constructed 
of wood and one of wire, was mounted in a wooden frame shown 


in Figure 5. [Ref. 4] The grid, placed 73 inches upstream of 
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the model pivot (Figure 6), creates nearly isotropic 
homogeneous’ turbulence. Table 1 summarizes measured 
turbulence intensities and estimated length scales determined 
by Roane. 

2. VLSAM Model and Mount 

The VLSAM model is a cruciform tail-control missile 
with low aspect ratio wings. The nose, wings, strakes, and 
tails are detachable from a hollow cylindrical body which can 
be rotated in 45 degree increments. Figure 7 portrays the 
VLSAM specifications. [Ref. 2] The model was rigidly mounted 
in the test section using a sting mount attached to a rotating 
arm which allows the angle of attack to be changed (Figure 6). 

3. Traverser and Five-Hole Pressure Probe 

The Velmex 8300 2 =D traverser used three 
microcomputer-controlled stepping motors mounted on top the 
test section to allow for manual or computer operated movement 
of an attached five-hole pressure probe (Figure 8). Traverser 
movement was controlled using software (PPROBE) and a 
microcomputer. 

Pressure measurements at each data point were acquired 
using a nulling five-hole pressure probe attached to the 
traverser (Figure 9). Pressure sample averages for each port 
were recorded at each data point and used to determine 
velocity, pitch angle, total pressure, total pressure 


coefficient, static pressure, and static pressure coefficient. 
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Figure 6. Planview of VLSAM Model with Pressure Probe and 
Grid in the Test Section of the Wind Tunnel (not drawn to 
scale) [Ref. 2] 


Intensity Length Scale Turbulence/ Dynamic Pressure 
(percent) (in.) Nfodel Dia. (Ih/ft?) 
Bye 1.05 eo 
0.89 | 1.88 
Three 0.62 16.38 
| O.15 15.01 
None . 15.85 





Table 1. Grid Turbulence Parameters [Ref. 5] 
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Figure 7. Specifications of VLSAM Model (Ref. 2] 
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Figure 8. 3-D Traversing Assembly 
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Figure 9. Five-Hole Pressure Probe [Ref. 22] 
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Probe port Pl measured stagnation pressure. Ports P2 and P3 
were static ports connected to a portable manometer. Flow yaw 
angle was read directly by rotating the probe until P2 and P3 
were equal. P4 and P5 were used with calibration curves to 
determine pitch angle. Velocities were found using the static 
and stagnation pressures as described previously. [Ref. 22] 
4. Data Acquisition 

A 48-port scanivalve attached to the five-hole probe 
allowed the pressures to be measured and converted to voltages 
by one transducer. A Hewlett-Packard data acquisition system 
consisting of a relay multiplexer, digital multimeter, relay 
actuator, and software enabled system control from a 


microcomputer. [{Refs. 6,7,24] 


C. EXPERIMENTAL CONDITIONS 

Test conditions used by Roane [Ref. 4], Lung [Ref. 5], 
Viniotis [Ref. 6], and Johnson [Ref. 7] were repeated to allow 
for data comparison. 

(1) Turbulence grid three was used to generate freestream 
turbulence. Rabang found that grid three generated vortex- 
length-scale turbulence which had the largest side force 
effect on the missile. [Ref. 2] 

(2) A Reynolds number of 1.1*10? was obtained using 
reference pressure differences of 7.2 cm H,O for no grid and 


10.0 cm H,0 for grid three. These settings were determined 


a2 


by Roane [Ref. 3] to insure comparable test section velocities 
for each tunnel configuration. These corrections are 
necessary because the pressure difference measured by the 
manometer includes an error resulting from disturbances caused 
by the screen. The calibration constants determined by Roane 
were used to calculate the actual dynamic pressure. 
(3) Three body configurations were used. (Figure 10) 

- Body A: wings and tails at O degree roll angle 

- Body B: no wings or tails 

- Body C: wings and tails at 45 degree roll angle 
Body B was used for pressure measurements and all 


configurations were used for flow visualization. 





Figure 10. VLSAM Model Body Configurations [Ref. 2] 


(4) The nose roll angle was set at position eight. 
Rabang [{Ref. 2] found that the missile model's nose roll angle 


affected the side force coefficient, as suggested by previous 
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investigators. He determined that position eight produced the 
largest side force. 

(Gs) Missile angle of attack was set at 50 degrees, 
determined by Rabang to be the angle of attack for the maximum 
Side force. [Ref. 2] 

(6) The data plane was at eleven body diameters from the 
nose corresponding to that used by Johnson. [Ref. 7] 

(7) Wind tunnel temperature was not allowed to increase 
by more than 20 degrees during a data run. When the 
temperature reached the limit the run was stopped until the 


temperature cooled, at which point the testing resumed. 


D. SOFTWARE 

PPROBE (Appendix A) iS a program written by Kindelspire 
(Ref. 25] and modified by Lung [Ref. 5] to control traverees 
movement. The operator may choose manual or computer control, 
the size of plane to be measured, and the step size for each 
movement. 

CALP (Appendix B) iS a program correlating transducer 
output voltage and pressure using a calibration manometer. 
The program was executed prior to and upon completion of the 
data collection run. CALP's output was converted to a slope- 
intercept equation using linear regression on a hand-held 
calculator. The two equations were averaged and the results 


used in CONVERT to calibrate the scanivalve transducer output. 
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CONVERT (Appendix C) is a program written by Lung [Ref. 


5] which reads the PPROBE data file and computes total 
velocity, x and y components of velocity, pitch angle, total 
pressure, total pressure coefficient, static pressure, and 


static pressure coefficient for each data point. Velocities, 
pitch angle, and pressures were calculated using calibration 


curves provided by the probe manufacturer. Pressure 


coefficients were calculated using ambient room temperature 


and tunnel calibration factar (K), and were non- 
dimensionalized by the dynamic pressure. Equations (3) and 


(4) were used to calculate the static and total pressure 


coefficients, respectively. 


C..= (P,i-P,) /Q (3) 
Coes (Prim Py) JQ (4) 
Where: 

Sas Static pressure coefficient 

C,, = Total pressure coefficient 

Q = Freestream dynamic pressure 

Pe > Hreestream static pressure 

Pee! Leecebean eocal) pressure 

Pe eecal Static pressure 

P,, = Local total pressure 


EL 


VORTIC is a program written by Johnson [Ref. 7] to compute 


the dimensionless vorticity using the velocity output from 


CONVERT using equation (5). (Appendix D) 


2D 


V= (0 (v/U0) 7o(x/ a) } =o (a7 U7 ae (5) 
Where: 
V"= Vortieity 


y component of the velocity 


< 
tl 


streamwise coordinate distance 


Ma 
tl 


u = X component of the velocity 


y transverse coordinate distance 
TECPLOT is a commercial graphics program used to generate 
the velocity vector, pressure coefficient contour, and 


VOrtlcitcy conteur pleter 


E. EXPERIMENTAL PROCEDURE 


The following experimental procedure was used to collect 


data. 

(1) The data acquisition equipment and wind tunnel were 
energized. Ambient tunnel temperature and pressure were 
recorded. The wind tunnel speed was monitored using the 


appropriate manual manometer value corresponding to the tunnel 
configuration previously discussed. 

(2) Pressure calibration was performed and the results 
recorded. 

(3) PPROBE was executed. The probe was positioned to the 
start. point using manual “ecencrol Computer control was 
selected followed by grid size, increment distance and input 


file name. The probe was physically rotated until the 
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Peterence between FP, and FP, wWaS approximately zero as 
indicated by the portable digital manometer. The yaw angle 
was read from the probe and entered on the computer. Each of 
the five pressure ports was automatically sampled ten times, 
averaged, and displayed on the computer screen. The values 
displayed were in units of volts and subsequently converted 
Eeepressure Units. The P,; measurement had to be positive or 
a very small negative (> -0.5 volts) value to ensure the probe 
memearcing the Correct direction. Ime difference between P, 
ema FF, Values had to be less than 0.1 volts, which was 
approximately a 6.4% error of the dynamic pressure for the no- 
grid configuration and an 8.6% error of the dynamic pressure 
meeche grid #3 configuration (the difference between P, and 
P, would ideally be zero). If either requirement was not 
satisfied, the probe was adjusted and the point resampled. 
At the completion of each column of data points PPROBE would 
automatically store the x-y position and corresponding average 
values and move the traverser to the next column. This 
procedure was repeated for the remaining data plane. 

(4) CALP program was executed, the results were recorded, 
and the calibration curve was calculated as previously 
described. 

(5) CONVERT and VORTIC programs were executed to obtain 


necessary data for plotting results. 


Za 


(6) Output files from CONVERT and VORTIC were converted 
to input file format for TECPLOT. Velocity vectors, pressure 


coefficient, and vorticity contours were then plotted. 


F. FLOW VISUALIZATION EQUIPMENT 

The flow visualization system was designed for the NPS 
low-speed wind tunnel by Chlebanowski [Ref. 26] and modified 
by Sommers [Ref. 27]. The system is depicted in Figure 11. 

Smoke generated external to the tunnel by a portable 
fog/smoke machine was piped to a 10-gallon plenum chamber 
which served to dampen surges of smoke from the machine and 
allow the smoke to cool slightly. The smoke exited the plenum 
chamber and entered the contraction section of the tunnel via 
flexible hoses and copper tubing. The smoke then entered the 
freestream tunnel air three feet ahead of the missile model. 

A laser sheet was used to illuminate a thin plane on the 
leeward side of the missile, allowing the flow phenomenon to 
be observed and recorded. The laser source was a class four 
Spectra-Physics, five watt, Argon-ion laser. The laser 
emitted a 1.25 mm blue-green beam through a Newport 
Corporation F-LFI laser fiber to a fiber optic handset, beam 
expander, and plano-cylindrical lens. The handset, expander, 
and lens were mounted to a manual traverse mechanism above the 
test section. The traverse allowed the laser sheet to be 


moved to any position along the missile. 
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Figure 11. Flow Visualization System [Ref. 27] 


The flow visualization runs were video-taped using a 
color video camera, 3/4-inch VCR, and color television 
monitor. Still photographs were taken of selected tape frames 
by the NPS Educational Media Department directly from the 


Monitor screen. 


III. RESULTS AND DISCUSSION 


The following sections discuss the velocity vector plots, 
total pressure coefficient contour plots, vorticity contour 
plots, and flow visualization photographs for the missile 
model configuration B (body-only) at eleven diameters. The 
results produced without using a turbulence grid are labeled 
"OB" and with turbulence grid three are labeled "3B", 

The plots generated from Lung's data are at a survey grid 
plane six diameters (10.5 inches) aft from the missile nose. 
Lung also used the body-only configuration, but used grid one 
for turbulence generation. 

The plots view the missile as it was mounted in the tunnel 
when observed in the direction of the tunnel freestream. All 
measurements are in inches (indicted by "), followed by the 
equivalent missile model base diameters (1.75 inches and 
indicated by d), and will be referenced from the leeward side 
of the missile body and the centerline. The terms "above" and 
"below" are relative to the missile centerline. "Away" 


indicates the distance from the center of the missile body. 


A. CONFIGURATION OB 
The velocity vector plot (Figure 12) indicates the 


location of three vortices and shows the location of a saddle 
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Figure 12. Velocity Vectors-Configuration OB 
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point. The two larger vortices are asymmetric with the flow 
in opposite directions. The upper one is located about 1" 
(0.57a) above and 2 1/8" (1.21d) away. The lower outside one 
is located 11/2" (0.86d) below and 6 1/4" (3.57d) away. The 
lower inner one is located 7/16" (0.25d) below and 2 1/4" 
(1.29d) away. The saddle point is located 5/8" (0.36da) below 
and. 3 1/8” (1.79a)-avayz 

The total pressure coefficient plot (Figure 13) shows one 
large center of pressure loss located on the centerline and 
4 1/4" (2.34d) away. The center does not coincide with any 
core from the velocity vector plot. The contours change from 
C,, level five to level two (-0.29 to -0.82) in 3/4" in the 
upper half and in 1" in the lower half. This larger gradient 
in the upper half indicates a greater vortex strength. 

The vorticity plot (Figure 14) shows three centers. The 
largest upper one is 5/16" (0.18d) above and 5 1/4" (3.00d) 
away. The smaller upper one, which is weaker, is 5/16" 
(0.18d) above and 4 1/2" (2.57d) away. The lower one is 1 
3/16" (0.79d) below and 4 7/8" (2.79d) away. The Zero vale 
alternates above and below the missile centerline with no 
apparent pattern. The upper vorticity cores have positive 
VOEE TC ikea The lower one has negative vorticity and is 
smaller in magnitude than either of the upper two. 

The vorticity plot from Lung (Figure 15) shows two 


centers. The upper one core is 3/4" (0..43da) above and 2 7/8" 
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Figure 13. Total Pressure Coefficient-Configuration OB 
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Figure 14. Vorticity-Configuration OB 
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Figure 15. Vorticity-Configuration OB (six diameters) 
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(1.64d) away. The lower one is 5/8" below (0.36d) "and 35a) 
(1.86d) away. The upper core has the largest magnitude and 
the greatest gradient. There also appear to be two centers 
at the inner edge of the survey grid, one above and one below 


thescentverlines 


B. CONFIGURATION 3B 

The velocity vector plot (Figure 16) indicates the 
location of three vortices and a saddle point, similar to the 
OB vector plot. The upper one is located 7/8" (0.50d) above 
and 2 1/8" (1.21d) away. The larger lower one is located 1 
1/2" (0.86d) below and 6 3/8" (3.64d) away. The smaller lower 
one is located 5/8" (0.39d) below and 3 1/8" (1.79d) away. 
The saddle point is located 5/8" (0.36d) below and 3 1/8" 
(1.79d) away. The two larger vortices are asymmetric with 
opposite flow direction. 

The total pressure coefficient plot (Figure 17) shows one 
large center of loss located 5/16" (0.18d) above and 4 5/8" 
(2.64d) away. The center does not coincide with any core from 
the velocity vector plot. The contours change from C,, level 
six to level two (-0.28 to -0.86) in 3/4" (0.43d) in the upper 
half and in 1 1/4" (0.71d) in the lower half. This larger 
gradient change indicates a greater vortex strength in the 


upper half. 
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Figure 16. Velocity Vectors-Configuration 3B 
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Figure 17. Total Pressure Coefficient-Configuration 3B 
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The vorticity plot (Figure 18) shows two centers. The 
upper one is 1/2" (0.29d) above and 5 1/2" (3.14d) away. The 
Mewer core 1s 1 1/78" (0.64d) and 5 1/4" (3.00d) away. The zero 
value is roughly along the missile centerline. The upper core 
has positive vorticity. The lower one has negative vorticity 
and is smaller in magnitude than the upper core. 

The vorticity plot from Lung (Figure 19) represents 
turbulence added by use of grid #1. The plot shows one center 
1/2" (0.29da) above and 3" (1.71d) away. A second center is 
indicated in the lower half. The upper core has a much higher 


magnitude and gradient than the lower one. 


C. COMPARISON BETWEEN NO TURBULENCE AND TURBULENCE 

The addition of turbulence on the scale of the vortices 
did not affect the overall shapes of the plots. The effect 
on the velocity vector plot was to decrease the magnitudes, 
but the directions and core positions were approximately the 
same. The addition of turbulence did not change the 
magnitudes of the total pressure coefficient contours, but the 
gradient was much tighter and the center of the core moved up 
pero” (0.18d). The vorticity strength was less and the 
gradients were greatly reduced with turbulence. The vorticity 
plots from Lung showed similar changes between the no 


mrroulence and turbulence conditions at six diameters. 
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Vorticity-Configuration 3B 
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Figure 19. Vorticity-Configuration OB (six diameters) 


41 


These results indicate that the added turbulence smoothes 
and steadies the flow, and results in a decrease in vorticity. 
This decrease in vorticity for added turbulence corresponds 
to Rabang's findings that the side force decreases with added 


turbulence for the body-only configuration [Ref. 2]. 


D. COMPARISON BETWEEN BODY-ONLY CONFIGURATIONS 

Comparison of the vorticity plots for the body-only 
configuration at six diameters and eleven diameters for the 
no turbulence condition shows a higher vorticity strength at 
eleven diameters, but the gradients are approximately equal. 
The velocity vector plots (Figures 12,20) show that the 
positions of the vortex cores are farther away from the 
missile body at eleven diameters although the relative 
positioning is about the same. The total pressure coefficient 
plots (Figures 13,21) show two tight cores and a larger 
gradient at six diameters compared to one large core at eleven 
diameters. The second core iS expected to be outside the 
survey grid. 

For the turbulence condition, however, the vortex strength 
1s much higher and the gradient is larger at six diameters. 
The velocity vector and total pressure coefficient plots show 
the same relative characteristics as the no turbulence 
Condesa ton: 

These results indicate that the added turbulence reduces 


the vorticity at both six diameters and eleven diameters, but 
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Figure 20. Velocity Vectors-Configuration OB (six diameters) 
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Figure 21. 


Total Pressure Coefficient-Configuration OB (six 
diameters) 
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the magnitude of the vorticity reduction is greater at eleven 


diameters. 


E. COMPARISON BETWEEN BODY-ONLY AND WINGED CONFIGURATIONS 
Comparison of the body-only ("0B" and "3B") configuration 
Geo the "x" COC Sandu oe jo and) +" COAT Aand "“3A™) 
configurations [{Ref. 7] at eleven diameters shows the body- 
only configuration to have a strong resemblance to the "x" 
configuration and little resemblance to the "+" configuration 
for both the no turbulence and turbulence conditions. 
Comparisons of the velocity vector plots between "0C" and "3c" 
(Figures 22,23) and "OB" and "3B" (Figures 12,16) show the 
Same positioning and flow directions of the vorticity cores 
and the saddle points. Comparison of the total pressure 
coefficient plots (Figures 13,17,24,25) show comparable 
coefficient magnitudes, but the body-only configurations have 
a much larger gradient. Comparison of the vorticity plots 
with the "x" configuration shows a higher vorticity for the 
body-only configuration with no turbulence (Figures 14,26) and 
comparable vorticity magnitudes for the turbulence condition 
Sergures 18,27). Comparisons with the iA eoOnrlqurat lems 
(Figures 28,29) show that the body-only configuration 
Magnitude is greater with no turbulence but is less for the 


mirtbulence condition. 
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Figure 22. Velocity Vectors-Configuration OC 
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Figure 23. Velocity Vectors-Configuration 3C 
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Figure 24. Total Pressure Coefficient-Configuration OC 
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Figure 25. Total Pressure Coefficient-Configuration 3C 
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Figure 26. Vorticity-Configuration 0C 
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Figure 27. Vorticity-Configuration 3C 
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Figure 28. Vorticity-Configuration OA 
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Figure 29. Vorticity-Configuration 3A 
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These results indicate that the addition of wings, 
strakes, and tails in the "x" configuration tends to reduce 
the magnitude of the vorticity with no turbulence. Again, as 
in the body-only configuration, the addition of turbulence of 
the nose-generated vortex scale appears to steady the flow and 
results in a decrease in vorticity for the "x" configuration. 
For the "+" configuration the turbulence has less of an effect 
on the vorticity magnitude, indicating a dependence on roll 
angle also. In addition, the amount of decrease in vorticity 


is larger for the body-only configuration. 


F. COMPARISON WITH FLOW VISUALIZATION 

Comparison of the body-only configuration velocity vector 
plots at six and eleven body diameters with flow visualization 
photographs (Figures 30,31) show that the velocity vector plot 
cores correspond to the asymmetric vortices formed on the lee 
Side of the missile. Figure 32 is a sketch illustrating the 
vortex core movement away from the missile as the distance 
from the missile nose is increased from six to eleven body 


diameters. 
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Figure 30. Flow Visualization (six diameters) 
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Figure 31. Flow Visualization (eleven diameters) 
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Figure 32. Sketch of Vortex Core Positions 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

The flowfield about a vertically-launched surface-to-air 
missile (VLSAM) model at high angle of attack was studied. 
The body-only configuration at eleven body diameters was 
tested and compared to the body-only configuration at six body 
diameters and with the winged configurations at eleven body 
diameters. The angle of attack was set at 50 degrees with a 
Reynolds number of 1.1*10? for all data runs. All tests were 
conducted with both an ambient wind tunnel condition (no grid) 
and a condition with grid-generated turbulence on the scale 
of the nose-generated asymmetric vortices. The following 


conclusions were reached: 


(1) Addition oof turbulence for the body-only 
configuration tended to smooth and steady the flow. The 
turbulence reduced the vorticity level, but did not 


Significantly change the positions of the vortices. 

(2) The addition of turbulence had a greater effect on 
vorticity reduction at eleven body diameters than at six body 
diameters for the body-only configuration, based on the 
vorticity plots (Figures 14) 153ieio 2 

(3) The body-only plots at eleven body diameters are very 
Similar to the "x" configuration at the same point, ~bue 
different from the "+" configuration. Without turbulence, the 


vorticity levels for the body-only case are much larger than 
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fon the "x" and the "+" configurations. With turbulence, 
however, the comparisons show that the vorticity levels were 
not only reduced but were approximately equal between the 
body-only and "x" configurations. However, compared to the 
"+" configuration, the vorticity levels were lower. 
Apparently the effect of added turbulence.on vorticity is much 
greater for the body-only configuration than the winged 
configurations. The winged cases are also complicated by a 
dependence on roll angle. Rabang noted the same results on 
the side force magnitudes. 

(4) The addition of wings had a significant effect on the 
reduction of the vorticity level without added turbulence. 
With turbulence, the addition of wings had a negligible 
effect on the vorticity level. Again, Rabang noted the same 
results on the side force magnitude. Apparently, the effects 
of added turbulence outweighs the effects of adding wings. 

(5) The flow visualization experiment verified that the 
asymmetric vortices move away from the missile as distance is 
increased from the missile nose along the body axis, in a vary 
non-linear manner. The velocity vector plots tended to match 
the flow visualization results for the location of the vortex 


cores. 
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Recommendations for future research are: 

(C) Conduct pressure measurements outside the present 
survey grid to determine characteristics of secondary vortices 
that may not have been located near the missile or farther 
away from the missile. 

(2) Attempt to decrease the vortex asymmetry by modifying 


the present missile nose geometry. 
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APPENDIX A. 


BEF SEG: CLEAR , &HFEOO: GOTO 
GOTO 1000 ' User program 
GOTO 900 ' Error handling 


PPROBE PROGRAM 


4 'Begin PCIB Program Shell 


I=&HFEOO ' Copyright Hewlett-Packard 1984,1985 


PCIB. DIRS=ENVIRONS ("PCIB") 
mo—PerB. DIRSt’\PCIBILC. BLD" 
BHOAD IS, I 
CALL I(PCIB.DIRS,1%,J%) :PCI 
IF J%=0 THEN GOTO 13 
PRINT "Unable to load."; 
PRINT " Chr DOtaey. dio 5 ‘ye! 
END 


B.SEG=1% 


BEF SEG=PCEBSSEG:0O.5—5:C.S=10:1.V=15 


fm.,C=20:L.P=25: LD. FILE=30 
GET .MEM=35:L.S=40: PANELS=4 
PCIB. ERRS=STRINGS (64,32) 


5:DEF.ERR=50 
PCIB.NAMES=STRINGS (16, 32) 


CALL DEF. ERR(PCIB. ERR, PCIB. ERRS, PCIB.NAMES, PCIB.GLBERR) 
PCIB. BASERR=255 


ON ERROR GOTO 3 
J=-1 
IS=PCIB.DIRS+"\PCIB.SYN" 
CALL 0.S(IS$) 
IF PCIB.ERR<>0 THEN ERROR 
I=0 

5 


PCIB. BASERR 


Cc A L L 


fee READ. REGISTER, READ. SELFID, DEFINE, INITIALIZE.SYSTEM) 
Poelr PCIB.ERR<>0 THEN ERROR PCIB. BASERR 


Z 


7 


C A L L 


fe), ENABLE .SYSTEM, DISABLE.SYSTEM, INITIALIZE, POWER.ON) 
met F PCIB.ERR<>O0 THEN ERROR PCIB. BASERR 
Semeenits: 1.V(L,MEASURE, OUTPUT,START, HALT) 
fee PCIB.ERR<>0 THEN ERROR PCIB. BASERR 


3 


if 


C A L L 


iol, ENABLE. INT. TRIGGER, DISABLE. INT.TRIGGER, ENABLE.OUTPUT,D 
PSABLE.OUTPUT) 


5 
3 
34 


IF PCIB.ERR<>0 THEN ERROR 
@20G L.V{L,CHECK. DONE,GET. 
IF PCIB.ERR<>0 THEN ERROR 
CALL I.V(1I,SET.MODE, WRITE. 
IF PCIB.ERR<>0 THEN ERROR 


PCIB.BASERR 
STATUS, SET. FUNCTION, SET. RANGE) 
PCIB.BASERR 

CAL, READ. CAL, STORE. CAL) 
PCIB.BASERR 


CALL I.V(I,DELAY,SAVE.SYSTEM,J,J) 


Pee CIB.ERR<>0 THEN ERROR 
I=1 


PeiB.BASERR 


CALL I.V(1,SET.GATETIME, SET. SAMPLES , SET. SLOPE, SET. SOURCE) 


onl 


41 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

42 CALL I.C(I,FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 
43 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

44 CALL I.C(I,INTERVAL, RATIO, TOTALIZE,R1OOMILLI) 
45 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

46 CALL 1.G(1,Ri, RO RO Omicron 

47 IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

48 CALL I.C(I,R10MEGA,R100MEGA, CHAN.A,CHAN.B) 
49 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

50 CALL I.C(I, POSITIVE, NEGATIVE, COMN, SEPARATE) 
51. IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

52 I=2 

53 I=3 

54 CALL I.V(I,ZERO.OHMS,SET.SPEED,J,d) 

55 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

56 CALL I.C(I,DCVOLTS ,ACVOLTS , OHMS , R200MILLTI) 
57 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

58 CALL I.C(1I,R2,R20, R200) Romer 

59 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

60 CALL I.C(I,R20KILO, R200KILO, R2MEGA, R20MEGA) 
61 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

62 CALL I.C(I,AUTOM,R2.5,R12.5,J) 

63 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 


64 I=4 
6 5 S A E ie 
I.V(I,SET. COMPLEMENT, SET. DRIVER, OUTPUT.NO.WAIT, ENABLE. HANDSH 
AKE) 

66 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

6 7 . A Ti L 


I.V(I,DISABLE. HANDSHAKE, SET. THRESHOLD, SET.START. BIT, SET. NUM. 
BITS) 

68 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

69 CALL IsV(1, SET. LOGIG oar sun 

70 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

71 CALL I.C(I, POSITIVE, NEGATIVE, TWOS , UNSIGNED) 

72 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

73 CALL I.C(I,0¢,TTLARon a 

74 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

75 CALL I.C(1,R2, Ragen 

76 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

77 CALL TsC(1-R6 Ri oR See 

78 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

79.CALL I.C(1I,R10, Riise 

80 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

81 CALL I.C(1I,R14,RiS seer 

82 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

83 I=6 

8 4 Cc A ve Ti 
I.V(1I,SET. FREQUENCY, SET.AMPLITUDE, SET. OFFSET, SET.SYMMETRY) 
85 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

86 CALL I.V(1I,SET.BURST.COUNT,J,J,J) 
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87 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

88 CALL I.C(1I,SINE, SQUARE, TRIANGLE, CONTINUOUS) 

89 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

So CALL I.C(I,GATED, BURST,J,d) 

91 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

92 I=7 

9 3 € A e if 
I.V(I,AUTOSCALE, CALIBRATE, SET. SENSITIVITY , SET. VERT. OFFSET) 
94 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

9 5 @ A I, ; 
I.V(1I,SET.COUPLING, SET. POLARITY, SET. SWEEPSPEED, SET. DELAY) 

96 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 


9 Zi Cc A id 1 
m7 (rt, SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET. TRIG. LEVEL, SET. TRIG 
. MODE) 

98 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

9 9 © A L i 


I.V(I,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO,GET.TIMEBASE. IN 
FO) 

100 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1 O 1 Cc A ia mi 
I.V(1I,GET.TRIG. INFO, CALC.WFVOLT, CALC.WFTIME,CALC.WF.STATS) 
102 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1 0 3 Cc A ie T 
I.V(I,CALC.RISETIME, CALC. FALLTIME, CALC. PERIOD, CALC. FREQUENCY) 
104 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 


1 0 5 C A L if 
I.V(I,CALC.PLUSWIDTH, CALC. MINUSWIDTH, CALC. OVERSHOOT, CALC. PRE 
SHOOT) 

106 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1 0 7 Ee A it it 


I.V(I, CALC. PK. TO. PK, SET. TIMEOUT, SCOPE. START , MEASURE. SINGLE. WF) 
108 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

109 CALL I.V(I,MEASURE.TWO.WF,J,J,J) 

110 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

111 CALL I.C(I,R10NANO, R100NANO, R1IMICRO, R10MICRO) 
112 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

113 CALL I.C(I,R100MICRO, RIMILLI,R1OMILLI, R1OOMILLT) 
114 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

115 CALL I.C(I,R1,R10,R20NANO, R200NANO) 

116 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

117 CALL I.C(I,R2MICRO, R20MICRO, R200MICRO, R2MILLT) 
118 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

jromcAT.T, £.C(l,R20MILLI, R2Z00MILLI, R2, R20) 

120 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

121 CALL I.C(I,R50NANO, R500NANO, R5MICRO, R50MICRO) 
122 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

123 CALL I.C(I,R500MICRO, RSMILLI, RSOMILLI, R500MILLT) 
124 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

125 CALL I.C(I,R5,R50,CHAN.A, CHAN.B) 

126 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 


63 


ee 
128 
29 
Ls0 
ibs gil 
1 
133 
134 
135 
136 
L37 
PS 
ay, 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
P50 
itis) Al 
E52 
153 
154 
i Pays: 
156 
oe, 
158 
Joye, 
160 
1G 
eZ 
16S 
164 
800 
801 
Sag 
900 
900 
902 


CALL I.C(I,EXTERNAL, POSITIVE, NEGATIVE, AC) 
IF PCIB.ERR<>0 THEN ERROR PCPS EASERR 
CALL I.C(I,DC, TRIGGERED, AUTO. PREG )AUTO. LEiaam 
IF PCIB.ERR<>0 THEN ERROR PCiepe ea en 
CALL I.C(I,X1,X10,STANDARD, AVERAGE) 
IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
I=8 

CALL I.V(I,OPEN. CHANNEL, CLOSE. CHANNEL,J,J) 
IF PCIB.ERR<>0 THEN ERROR PCI 2B SASSER: 
CALL C.S 

IF PCIB.ERR<>0 THEN ERROR PCIB. BASERR 
ITS=PCIB. DIRS+"\ reuse. eee 

GALL Le (ro 

IF PCIB.ERR<>0 THEN ERROR SPGfe. SASERE 
IS="DMM. 01": I=3 :J=0: K=0: L=1 

CALL DEFINE(DMM.01,1$,1,J,K,L) 

IF PCIB.ERR<>0 THEN EBRRORSPere eaoenk 
IS="Func.Gen.01":I1=6:J=0:K=1:L=1 

CALL DEFINE (Func. Gen .@1) 05,4) 0.,.-92) 
IF PCIB.ERR<>0 THEN ERROR P@re eee. 
IS="Scope.01": I=7:J=0: K=2: L=1 

CALL DEFINE (Scope.01) 1s) Dmkem) 

IF PCIB.ERR<>0 THEN ERROR P@lb eases n. 
TS="Counter.01": I=1:J3=0: K=3: L=1 

CALL DEFINE(Counter sol, lon) aiitae 

IF PCIB.ERR<>0 THEN ERROR PG@is, becoerre 
IS="Dig.In.01": 1I=4:J=0:K=4:L=1 

CALL DEFINE (Dig. Ins Ouse ce ese Ken 

IF PCIB.ERR<>0 THEN ERROR PCB SBA5ERE 
IS="Dig.Out.01": T=4:J=1: K=4: L=1 

CALL DEFINE(Dig sour. Ol 10> paren 

IF PCIB.ERR<>0 THEN ERROR PCIE BASERE 
IS="Relay.Act.01":I=8:J=0:K=5:L=1 
CALL DEFINE(Relay.Act.01,1$,1,J,K,L) 
IF PCIB.ERR<>0 THEN ERROR PCie sacar 
IS="Relay.Mux.01":I=2:J=0:K=6:L=1 
CALL DEFINE(Relay.Mux.01,iIS7f,d,, 5) 
IF PCIB.ERR<>0 THEN ERROR PCB EA eee 
IS=ENVIRONS (""PANELS") +"\ PANELS. EXE" 
CALL Lwsiers) 

GOTO 2 

IF ERR=PCIB.BASERR THEN GOTO 903 
PRINT "BASIC error #";ERR;" occurred in line ";ERL 
STOP 


903 TMPERR=PCIB.ERR: IF TMPERR=0 THEN TMPERR=PCIB.GLBERR 

904 PRINT "PC Instrument error #";TMPERR;" detected at line 
My ERs 

905 PRINT "Error: “;P@rPBsEer-— 

906 IF LEFTS$(PCIB.NAMES,1)<>CHR$(32) THEN PRINT "Instrument: 
">PCIB.NAMES 

907 STOP 
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908 COMMON PCIB.DIRS$,PCIB.SEG 
909 COMMON LD.FILE,GET.MEM, PANELS, DEF.ERR 


g 1 8) Ss O M M O N 
PCIB. BASERR, PCIB. ERR, PCIB. ERR$, PCIB.NAME$, PCIB.GLBERR 
2, 1 i Cc O M M O N 


READ. REGISTER, READ. SELFID, DEFINE, INITIALIZE.SYSTEM, ENABLE.SY 
STEM, DISABLE.SYSTEM, INITIALIZE, POWER.ON, MEASURE, OUTPUT, START 
Aaa pENABLE.INT.TRIGGER, DISABLE. INT.TRIGGER, ENABLE.OUTPUT,D 
moo be OUTPUT, CHECK. DONE,GET.STATUS 912 COMMON 
SET. FUNCTION, SET. RANGE, SET.MODE, WRITE. CAL, READ. CAL, STORE. CAL 
memiaAyY ,oAVE.oYSLEM,SET.GATETIME,SET.SAMPLES,SET.SLOPE,SET.SO 
Mreh ZERO.OHMS ,SET.SPEED, SET.COMPLEMENT, SET. DRIVER, OUTPUT.NO 
.WAIT, ENABLE.HANDSHAKE, DISABLE.HANDSHAKE 913 COMMON 
fae the oOnOnD SET. START. BIT, SERANUM. BITS ,SET.LOGIC.SENSE,SET 
Peewee y, oS) sAMPLITUDE ,SET.OFFSET,SET.SYMMETRY,SET.BURST.C 
eee TOSCALE, CALL BRATE,SEL.SENSITIVITY,SET.VERT.OFFSET,SET 
BeOUPLING, SET. POLARITY,SET.SWEEPSPEED 914 COMMON 
Peo ein y por TS TRIG. SOURCE,SET.TRIG.SLOPE,SET.TRIG.LEVEL,SET. 
eo, MODE, CEPTS SoINGCLE.WF,GET.IWO.WF,GET.VERT.INFO, GET.TIMEBAS 
Been PO, GET. TRIG. INFO, CALC.WFVOLT, CALC.WFTIME,CALC.WF.STATS,C 
ALC.RISETIME, CALC. FALLTIME, CALC. PERIOD Ons COMMON 
CALC. FREQUENCY, CALC. PLUSWIDTH, CALC.MINUSWIDTH, CALC.OVERSHOOT 
perv oc. PRESHOOT, CALC. PK.TO.PK,SET.TIMEOUT, SCOPE.START ,MEASURE 
-SINGLE.WF,MEASURE.TWO.WF, OPEN.CHANNEL, CLOSE. CHANNEL 2) ILS 
c O M M O N 
FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER, INTERVAL, RATIO, TOTALIZGE, 
R1OOMILLI,R1,R10,R100,R1KILO, R1OMEGA, R100MEGA, CHAN.A,CHAN.B, 
POSITIVE, NEGATIVE, COMN, SEPARATE, DCVOLTS , ACVOLTS , OHMS , RAOOMIL 
mei, RZ0,RZ00 ,R2KILO, R2Z0KILO, RZ00KILO pe ial ll 7 COMMON 
R2MEGA, R20MEGA,AUTOM,R2.5,R12.5, POSITIVE, NEGATIVE, TWOS, UNSIG 
NED,OC,TTL,RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R14 
,R15,R16,SINE, SQUARE, TRIANGLE, CONTINUOUS, GATED, BURST, R1O0NANO 
,R1OONANO, RIMICRO,R1OMICRO, R1OOMICRO 918 COMMON 
R1IMILLI,R1OMILLI,R1OOMILLI,R1,R10,R20NANO, R200NANO,R2MICRO,R 
Pert eCRO, RZ00MICRO, R2MILLI, REOMILLI ,R200MILLI,R2,R20,R5ONANO, 
RSOONANO, R5MICRO, RS5OMICRO, R500MICRO, RSMILLI,RSOMILLI,R500OMIL 
meek, R50, CHAN.A, CHAN.B, EXTERNAL, POSITIVE VALS. COMMON 
NEGATIVE,AC, DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL, X1,X10,STANDAR 
D, AVERAGE 

9 2 0 c O M M O N 
fol, Func.Gen.01,Scope.01,Counter.01,Dig.In.01,Dig.Out.01, 
Relay.Act.01,Relay.Mux.01 999 'End PCIB Program Shell 

1000 REM This step initialzes the HP system 


1@10 CLS 
mezQ OPTION BASE 1 
1 O 3 O D ag M 


P(5),PA(50,5),PP(50,5) ,XPT(50) , YPT(50) ,X(50) , Y(50) , YAW(50) 
1040 REM 

1050 CALL INITIALIZE.SYSTEM(PGMSHEL. HPC) 

1060 REM 
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1070 REM SET FUNCTIONON THE '‘'DMM'" , "RELAY MUX 222 

ACTUATOR! 

1080 REM 

1090 CALL SET.FUNCTION(DMM.01,DCVOLTS) 

1100 CALL SET.RANGE (DMM. 01,AUTOM) 

1110 CALL DISABLE.INT.TRIGGER(DMM.01) 

1120 CALL ENABLE.OUTPUT (RELAY.MUX.01) 

1130 CALL ENABLE. OUTPUT (RELAY.ACT.01) 

1140 REM **®¥*¥k®¥k¥KKKKKKK PROGRAM TRAVERSE ®®RRRKKRKKKKKKKKEKKESR 

PT 50 7REM 

1160 REM OPEN THE COM PORT AND INITIALIZE THE MOTOR 

SETTINGS 

1170 OPEN "coml:1200,n,8,1,s,¢s, 4s, ca ome 

1180 REM SET MOTOR DEFAULT VALUES 

1190 DATA 2000,2000,2000,2,2,2,0.000125, 0.000125, 0.000 

1200 READ V1,V2,V3,R1,R2, R37 eee ee 

1210 REM DEFINE CHARACTERS FOR DATA REDUCTION ALGORITHM 

1220 RN2S=""RENAME A:RAW.DAT " 

1230 HEAD1S = " # X i Pl P2 P3 

P4 P5 YAW !" 

1240 FORMATS= "## ##.## ##.#4 ###.#7#2% F###.F## F##. 477 

HHH HHH FFF. HHH FFF. FE" 

1250 PRINT 

1260 PRINT "kk KK KKK KKK KKK KKK KKK KKK KKK KEKE EM 

1270 PRINT “** USER MUST SELECT "CAPS LOCK' FUNCTION 2. 

1280 PRINT "Xk KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KKKKKKEEN 

1290 REM DISPLAY MOTOR DEFAULT SETTINGS 

1300 PRINT *" KEKKEKEKEKKKKKKKEKKEKKEKEKKKEKKEKKKKEKEKKKEEEMNM 

1310 PRINT " INITIALIZED VALUES FOR ALL MOTOR 

SELLING S 3. 

1320, PRINE VELOCITY = 1000 STEPS /siee 

1330 "PRINTS" RAMP(MOTOR ACCELERATION) = 2 (6000 

STEPS/SEC*2)" 

1340 PRINT " DEFAULT INCREMENTAL UNITS ARE 

TNeHEsS" 

1350 PRINT " KKK KK KKK KKK KKK KK KKK KKK KK KKK EEK 

1360 PRINT 

1370 PRINT "NOTE!! USE MANUAL CONTROL TO INITIALIZE PROBE 

POSITION BEFORE" 

1380 PRINT " SELECTING COMPUTER CONTROLLED MOVEMENT. 
te 

1390 PRINT 

1400 INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER “MAN? 

Or wep) CONS 

1410 IF CONS="CP" THEN 3490 

1420 REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY eae 

ACCELERATION RAMP 

1430 PRINT 

1440 PRINT 

1450 PRINT " DO YOU WANT TO CHANGE THE VELOCITY Sa] 


ACCELERATION RAMP" 
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1460 
1470 
1480 
109 Vo 
1490 
1500 


PRINT " DEFAULT SETTINGS? (Y or N)" 
PRINT 
PRINT "IF 'NO', THIS PROGRAM WILL THEN LET YOU DEFINE 


Pave Dro tANGH YoU WANT FO MOVE (IN INCHES). IF 'YES'," 
PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS 


FOR ANY MOTOR." 


> 10 
1520 
3 0 
1540 


PRINT 
PRINT 
PRINT 
INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? 


feeor N)";DS 


TS 50 
1560 
e570 
5c 0 
290 
1600 
1610 


foo" Ve TrieN 1590 
IF DS="N" THEN 2220 


REM 

REM *kxk*k OQPERATOR SELECTED MOTOR VARIABLES ***xx* 
PRINT 

PRINT 

Peer "“WwHilCH DEFAULT VALUE? (ENTER ‘1*FOR®VELOC OR '2' 


FOR ACCEL RAMP) ";L 


1620 
1630 
N)" 

1640 
1650 
1660 


Gh GOTO 1690,1930 
PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? (Y OR 


INPUT V$ 
IF VS="Y" THEN 1690 
PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION 


RAMP? (Y¥ or N)" 


1670 IF RS="Y" THEN 1990 

1680 IF RS="N" THEN 1450 

1690 PRINT 

1700 PRINT 

im_enoe INPUT “WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2, 
res ) '' ? 

1720 ON J GOTO 1730,1830,1880 

1730 PRINT 

1740 PRINT 

meeo INPUT “ENTER DESIRED VELOCITY OF MOTOR #1";V1 
1760 PRINT 

1770 PRINT 

1780 PRINT 

1790 PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? 
(Y OR N)" 

1800 INPUT VS 

1810 IF VS="Y" THEN 1690 

1820 IF VS="N" THEN 1430 

1830 PRINT 

1840 PRINT 

1850 INPUT "ENTER DESIRED VELOCITY OF MOTOR 2";V2 
1860 PRINT 

1870 GOTO 1780 

1880 PRINT 
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1890 
1900 
1910 
19210 
1936 
1940 
P20 
2, 
1960 
13776 
1980 
136 
2000 
2010 
2020 


PRINT 
INPUT 
PRINT 
GOTO 1780 

PRINT 

PRINe 

INPUT "WHICH MOTOR ACCEL RAMP DO YOU WANT TO CHANGE? (1, 


"ENTER DESIRED VELOCITY OF SOTO. 73 3. = 


Or yoy; K 


ON K GOTO 970, 206002 1Z0 
PRINT 
PRINT 
ENPUT 
PRINT 
PRINT 
PRINT 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #1";R1 


"DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER 


MOTOR? “{Y Grupa 


2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2456 
2160 
2270 
2180 
2990 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
227.0 
2 


INPUT RMS 
IF RMS="Y oe THEN S30 
IF RMS="N" THEN 1450 
PRINT 
PRINT 
INPUT 
Peat 
PRINT 
GOTO 2000 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
GOTO 2000 
REM 
REM DEFINE DISTANCE “LO °MOy E16 fer 
PRE 
PRINT 
PRINT 
REM INITIALIZE MOTOR INCREMENTS TO ZERO 
I1=0 
I2=0 
I3=0 
PRINT 

Z 8 O Jee R I N 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #2";R2 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #3";R3 


T tt 


KKEKEKKKKKKKEKKKKKEKKKEKKEKKK KKK KEKE KKKEKKKEKKKKKKKKKKKKKKKKEKKKEEEEN 


2290 PRINT " *x DEFINE WHICH MOTOR YOU WANT TO MOVE 
kM 
2300 PRINT " ** 
KN 
2310 PRINT " ** NOTE!!! A POSITIVE ('+') INCREMENT TO 
A MOTOR *x" 
2320 PRINT " ** MOVES TRAVERSER AWAY FROM THAT PARTICULAR 
MOTOR aM 
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2330 


2340 


THE FLOW 


2350 


WINDOW 


2360 


DOWNWARD 


2 


PRINT " ** 
KN 


PRINT " ** 
kaM 


PRINT " ** 
aN 


Pree ™ 


-- MOTOR #1 MOVES THE PROBE UPSTREAM AGAINST 
-= MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS 


k* -— 
*& & Ut 


3 7 0 E R af N T E 


MOTOR #3 MOVES THE PROBE VERTICALLY 


KHKKKHKKHKKAHK AKER KKK 


2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
25 1.0 
2520 
2530 


PRINT 
PRINT 
PVPOr "WHEGH MOTOR DO YOU WANT TO MOVE? (1,2, 
ON, GOLTG 2420, 2680,2970 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 
PRINT " **¥*xxxkkkk (ENTER DISTANCE IN INCHES) *****eeex" 
ENEOUT Il 

PRINT 
PRINT" 
PRINT 
PRINT 
PRINT " 
SRINT * 
PRINT 


ore 3) 


KRAEEKKKEKEKKEKEKKEKKKEKEKEKKEKEKEKKEKEKKKEKEEEEEMN 


"SUMMARY OF OPERATOR INPUTS:" 
MOTOR #1 VELOCITY = ";V1 
ACCELERATION RAMP = ";R1 
INCREMENTAL DISTANCE = 


moll; INCHES" 


2540 
2509.0 
N)" 

2560 
2570 


PRINT! XK KK KKK KKK KKK KKK KERR REE 
PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 


PRINT 


Peay “ENTER “N' TO START MOTOR MOVEMENT. ENTER ‘'Y' TO 


RETURN" 


Z5o0 
2590 
2600 
2620 
2620 
2630 
2640 
2650 
2660 
267 0 
2680 
Zo90 
2700 
2710 
27 20 
2/30 
2740 
2750 


WetNte bo  VARTABLE SELECTION SUBROUTINE." 

INPUT V$ 

IF VS="Y" THEN 1430 

GOSUB 3410 

PRINT 

PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?" 
INPUT C$ 

DP ees="Y" THEN 2220 

PeeCo="N" THEN 3260 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #2?" 

PRINT '"" xxxkxkkkk*k (ENTER DISTANCE IN INCHES) **** ex exe" 
ENPOT l2 

PRINT 

PRINT 


REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
PRINT KX RK KKK KK KKK KKK KKK KEKE KEKE EEE 
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27 6OR PRINT 

2770 PRINT "SUMMARY OF OPERATOR INPUTS: " 

2730 ERENT ” MOTOR #2 VELOCE ae 

2/905 P Rint” ACCELERATION RAMP = ® Ez 

2800 PRINT " INCREMENTAL DISTANCE = 
“27; eenleGHES | 

2810 PRINT! Ke RK KEK KKKKKEKEKEKKKKKEKEKEKKKKEKEKKEEEKN 

202 Um ere 

Zoe 0) PRINT 

2840 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y oF 
N)" 

2650 PRINT 

2860 PRINT “ENTER 'N' TO START MOTOR MOVEMENT. ENTER) Yous 
RETURN" 

2870 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

2830 INPULVys 

2890 IF VS="YUS THEN 1436 

2900 GOSUB Ss cero 

2910 PRENG 

2920 PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?”® 


229 310 
2940 
29510 
2960 
27 0 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
300 
3060 
3070 
3080 
3090 
HeT3 
3100 
3 EO 
3120 
3130 
3140 
N)" 

3150 
JaGO 


INPUT “Gs 

TR CS="y" THEN 2220 

IF CS="N" THEN 8260 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #3?" 


PRINT " *¥xkkkkkkk (ENTER DISTANCE IN INCHES) ****xxxxx" 
TN PU Daas 

PROGNE 

PRINT 

REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 


PRINT!!! KERR KKK KEKKEKKEKEKKEKKEKKEKEKEKEKKEKEKKEKEKEEEEEN 


PRINT 
PRINT 
PR" 
PRINTS 
PRINT = 


"SUMMARY OF OPERATOR INPUTS:" 
MOTOR #3 VELOCITY = "72 
ACCELERATION RAMP = "';R3 
INCREMENTAL DISTANCE = 


7 SENGHE So 


PRINT 
PRINT 
PRINT 
PRINT 
PREND 


KKEKEKKKEKKKEKKEKKKEKKKKEKKEKKKKKEKKEKEEEEN 


"DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 


PRINT 


PRINT “ENTER 'N' TO START MOTOR MOVEMENT. ENTER Yea 


RETURN" 


3170 
SS 6 
3150 
3200 


PRINT "TO VARIABLE SELECTION SUBR¢ UI ee. 
INPUT V$ 

IF VS="Y" THEN 1430 

GOSUB 3410 


iLO 


3210» PRINE 
S82 Z0P PRINT 


o230 INPUT 


"DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT 


yaser I MS 


3240 IF MS="Y"" THEN 2210 

3250 PRINT 

3260 PRINT "DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR 
MOVEMENT?" 

3270 PRINT " KKKEKKKKKK NOTE! ! 1 kee RKKKKKK 

3280 PRINT " ALL PREVIOUS MOTOR INCREMENT INPUTS HAVE BEEN 
ZeROLZED." 

3290 PRINT "PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED 
MOVEMENT." 

3300 PRINT "x*xx*x (TIF 'NO!', THE PROGRAM WILL END). *****" 
3310 PRINT 

3320 INPUT "DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y 
omen)" :NS 

fee30 ITF NS="V" THEN 3500 

3340 PRINT 

S550 PRINT 

3360 PRINT 

3370 PRINT mt KEEN 

3380 PRINT " THE PROGRAM HAS ENDED." 

3390 PRINT it KEEN 

3400 END 

3410 REM xxkkKkKKK MOTOR MOVEMENT SUBROUTINE ****x*x**xxkx 
3420 PRINT #1, supe ee ;> PRI Ne we i 
eee "Cl es C2="-C23"23C3="°C3 

Por PRINT #12, "L1I="sT1is"2:Vl="sV1;":Rl=";R1;> 

Seo PRINT #1, “"sT2="sT2s"sV2="sV2:"3:R2="3R2 

Poo PRINT ¢1, "T3="7133"3V3="3V37"2:R3=";R3;":@e" 

3460 RETURN 

3470 REM KRREEKKEKKKEKEKKKEKKKKEKRKRKKKERRKKKRKRKRKRRKRKRKRKRKRRKR RRR 

3480 REM KKK KKK KKK KKK KKK KKK KKK KKK RRR KKK KERR EKER 

3490 PRINT 

3500 REM ******x** COMPUTER CONTROLLED MOVEMENT ****xxxx 

3510 PRINT 

3520 PRINT "THE PRESSURE DATA WILL BE WRITTEN TO FILES ON 
PeaveE ‘A’ " 

3530 PRINT 

3540 PRINT "YOU WILL BE ASKED TO INPUT FILE NAMES FOR THESE." 
3550 PRINT 

3560 INPUT "IS A FORMATTED DISK IN DRIVE 'A'? PRESS 'ENTER'! 
MOnecCONTINUE" ;DS 

3570 PRINT 

3580 PRINT 

3590 PRINT 

3600 PRINT " KKEKEKKEKKKEKKKEKREKEKKEKKKEKRKEKKEKREKKRKEKEKKEREREEM 
3610 PRINT " xx NOtE 112 aM 
3620 PRINT " xx COMPUTER CONTROLLED MOVEMENT xe 
3630 PRINT " *x* IS PROGRAMMED WITH A ke 


7h 


3640 
31650 
3660 
3670 
3680 
3690 


xx DEFAULTED NEGATIVE MOTOR INCREMENT *x"! 
xx (i.e. MOTOR #3 WILL MOVE UPWARD eM 


**k “BY ENTERING A’ (7) 3" @is fanieee eM 
KHKKKKKK KKK KKK EKER KEKE KERR EN 


Peer" 
PEENT ” 
PRINT?" 
Jerse Ie 
PRINT 

REM SET INITIAL MOVEMENT DISTANCE AND NUMBER OF DATA 


POUNTSe LO Z2ZERe 


3700 
3710 
Cae) 
3730 
3740 
39 S10 
3760 
3770 
3780 


YOU WANT TO MEASURE." 


3790 
3800 


MOVE. 


3810 
3820 
3050 
3840 
3650 


HT=0 
WD=0 

DIST=0 

XPT=0 

YPT=0 

N=0 

PRINT 

PRINT 

INPUT "WHAT IS THE DIMENSION ( X , 
;WD, HT 


Y ) (IN INCHES) THAT 
PRINT 

INPUT 
Do 
YPT=INEQHE” DE Say) Set 
XPT=INT(WD /DIST)+ 1 

N=XPT*YPT 

PRINT 
PRINT 


"WHAT IS THE STEP (IN INCHES) THAT YOU WANT TO 


"THERE ARE ";XPT;" * "“;yYPpT;" = ";N;" POINTS sian 


MEASURED " 


3860 
3870 
3880 
3890 
3200 
3 oaEO 
3920 


NAMES 


3330 
3940 
395 0 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 


PRINT 

INPUT "ARE THE NUMBER OF POINTS IS OK.(Y OR N)";C$ 

IF CS="N" THEN 3780 

CLS 

N=XPT 

IF (N < 1) OR (N > 99) GOTO 3780 

REM *** GENERATING STRING STRING SEGMENTS FOR DATA FILE 


REM ** STRING NUMBER "1" 
REM ** ENDING STRING NUMBER 


BS = MIDS$(STRS$(1), 
ES = MIDS(STRS(N), 
XS = "XXXXXxX" 

EXS = ".DAT" 

CLS 
PRINT 
PRINT 
PRINT (X$ + BS + EXS); " 
PRINT 
PRINT 
INPUT 


Phe) 2 
Bans 


LS 


"DATA FILES WILL BE INCREMENTED FROM:" 


To ": (XS + ES + ee 


PEN Tie DATA FILE NAME (6 CHARACTERS MAX -- NO 


EXTENSION) ";F2$ 


4040 
4050 
4060 
4070 


PRINT 

PRINT 

IF LEN(F2$) > 6 OR LEN(F2$) < 1 GOTO 4030 
els 
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4080 


PRINT N; "DATA FILES WILL BE GENERATED AND INCREMENTED 


AS FOLLOWS:" 


4090 
4100 
4110 
4120 
4130 
4140 
N) ee 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4302 
4304 
4306 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4405 


PRINT 
PRINT 
se sUIN ak 
PRINT 
PRINT 
Peer 
C$ 

IF cs _ Ne 
IF cs = uy 
GOTO 4140 
CLS 

PRINT 
PRINT 

REM SET INITIAL POSITION DATA 

%(1)=-DIST 

Y(1)=-DIST 

FOR IX=2 TO XPT+1 

X (IX) =0 

NEXT IX 

HOR Y=2 TO .YPT+1 

Y(JY)=0 

NEXT JY 

FOR I=1 TO XPT 

I1=0 

I2=0 

I3=0 

rORS3—1 TO YPT 

REM MOTOR CP-CONTROLLED MOTOR MOVEMENT 

I1=0 

I2=0 

I3=0 

REM EACH POINT TAKE 10 TIMES READINGS 

X(I+1)=X(1)+DIST 

XPT (J) =X (I+1) 

¥(J+1)=Y(J)+DIST 

YPT(J)=Y(J+1) 

iru " ADJUST THE WHEEL TO MAKE THE P2 =P3, INPUT THE YAW 


(F2S + BS + EXS); " To "; (F2$ + ES + EX$S) 


"ARE THE NUMBER OF POINTS AND FILE NAMES OK. (Y OR 


GOTO 3780 
GOTO 4180 


ANGLE" ; YAW (J) 


4408 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 


PRINT 

INPUT " PRESS 'ENTER! TO START THE MEASUREMENT" ; MOVES 
REM 

REM READ FIVE CHANNELS AND DISPLAY THE DATA 
REM 

STEPPER=4 

SWITCH = 3 

HOMER=8 

BEEAYI = .1 

DELAYS = 1 

REM SET THE S.V PORT TO #4 


eS 


4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
47-00 
4710 
a2 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4 


FOR l=]. Tom. 
THYME = TIMER 

CALL OUTPUT (RELAY.ACT.01,STEPPER) 
CHKTIME = TIMER 

IF CHKTIME < (THYME + DELAY1) GOTO 4540 
CALL OPEN. CHANNEL(RELAY.ACT.01,SWITCH) 
CLS 

NEXTM@EL 

PRINT 

PRINT " NOW IS POINT ";J 

REM START MEASURE FROM PORT 4 TO PORT 8 
FORT Id=—1) Ooms 

CALL OUTPUT (RELAY.ACT.01,STEPPER) 
CHKTIME = TIMER 

IF CHKTIME < (THYME + DELAY2) GOTO 4640 
REM EACH PORT SAMPLE 10 TIMES 

FOR I[il= On ho 

ROUT=1 

CALL OUTPUT (RELAY.MUX.01, ROUT) 

CALL MEASURE (DMM.01, VOLTS) 

PA (it, Jd) =—VOEES 

NEX®? 22 | 

CALL OPEN. CHANNEL(RELAY.ACT.01,SWITCH) 
IF JJ=5 THEN 4760 

NEXT JJ 

REM HOME THE S.V PORT TO #48 

CALL OUTPUT (RELAY .ACT 201 ,HOMER) 

CALL OPEN.CHANNEL(RELAY.ACT.01,HOMER) 
REM 

REM DISPLAY THE SAMPLE DATA 

REM 

PRINT HEAD1$ 

FOR IS= 1 TO 10 

8 4 0 P R if N TP U S I N G 


FORMATS ;IS,XPT(J) , YPT(J) , PA(IS,1) , PA(IS,2) , PA(IS,3) , PA(1S, 4) 
,PA(IS,5),YAW(J) 


4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
5000 


NEXT IS 

REM 

REM AVERAGE THE DATA 

REM 

FOR JA = 19f0 5 

TOTAL = 0 

FOR IA = 1 TO 10 

TOTAL = TOTAL + PA(IA,JA) 
NEXT IA 

AVERAGE = TOTAL /10 

P (JA) =AVERAGE 

NEXT JA 

PRINT 

PRINT "THE AVERAGES ARE: " 
PRINT HEAD1$ 
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5010 FOR JD=1 TO 5 

5020 PP(J,JD)=P(JD) 

5030 NEXT JD 

5 0 4 O Je R I N a U S ele N G 
FORMATS ;J,XPT(J),YPT(J),PP(J,1),PP(J,2),PP(J,3),PP(J,4),PP(d 
, 5) ,YAW(J) 

5045 PRINT 

5048 PRINT USING "THE NULLING ERROR IS +#.####";PP(J,3)- 
PP(J,2) 

5049 PRINT 

5050 PRINT "DO YOU WANT RE-MEASURE AGAIN (Y / N)" 

5060 PRINT 

5062 PRINT "IF '‘Y' WILL RE-SAMPLE AGAIN." 

5064 PRINT 

5070 INPUT "IF 'N' WILL MOVE THE TRAVERSER STEP UPWARD (WAIT 
feo eC )":Cs 

5075 PRINT 

5080 IF CS$="Y" THEN 4405 

5082 IF CS="N" THEN 5090 

5084 GO TO 5070 

5090 IF J=YPT THEN 5160 

5100 REM 

5110 REM MOVE THE TRAVERSER STEP UPWARD. 

5120 REM 

5130 I3=-DIST 

5140 GOSUB 3410 

B50 NEXT J 

5160 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

5170 OPEN "A:\RAW.DAT" FOR OUTPUT AS #2 

5180 PRINT #2 ,HEAD1$ 

5190 FOR ID=1 TO YPT 

pee .6|COOttiéiOd P R I N T # 2 OP > er re 


FORMATS ; ID, XPT(ID) , YPT(ID),PP(ID,1),PP(ID,2),PP(ID,3),PP(ID, 
4),PP(ID,5),YAW(ID) 5210 NEXT ID 


5220 
2230 
5240 
a2 00 
2)) 

5260 
a2 7/70 
5280 
3290 
5300 
5310 
Bo 20 
5 30 
5340 
5250 
5360 


CLOSE #2 

REM *** GENERATING INCREMENTED DATA FILE NAME 

IF (I > 10) OR (I = 10) THEN IS = MIDS$(STRS(I), 2) 

IF (I < 10) THEN IS = (MIDS(STRS(0), 2) + MIDS(STRS(I), 


pazee—= (F25 + 1S + EXS) 
PRINT 
PRINT " WRITING DATA FILE "; FI2$ 


DF2$S=RN2$+F1I2$ 

REM ** RENAME DATA FILE 

SUE EG | Dae 

REM 

PoieMOVvVE THE PRAVERSER TO THE NEXT SAMPLE POSITION 
REM 

PRINT 

ir I=XPT THEN 5430 


Wie. 


5370 INPUT "THEN PRESS 'ENTER' FOR NEXT COLUMN, SAMPIE@ 905 2a 
MOVES 


5320 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 


I2=-DIsST 

I3=HT 

GOSUB 3410 

NEXT I 

CLS 

PRINT "ALL MOVEMENTS COMPLETE" 

ee LIN Ag 

PRIN 

PRINT "YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER 


MOVEMENT (Y OR N)?" 


5480 
5490 


PRINT 
PRINT "IF 'Y', THE PROGRAM WILL TAKE YOU TO MANUAL 


CONTROEVSUBRCUTINE. | 


5500 
5D 10 
SoZ 0 
55 a0 
5540 


PRINT “IF 'N', THE PROGRAM WILL END." 


PRINT 

INPUT "ANOTHER MOVEMENT";RS 
IF RS = "Y" THEN A370 

IF R$ = "N" THEN 3370 
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APPENDIX B. CALP PROGRAM 


DEF SEG:CLEAR ,&HFEO0:GOTO 4 'Begin PCIB Program Shell 
GOTO 1000 ' User program 
GOTO 900 ' Error handling 
I=&HFEOO ' Copyright Hewlett-Packard 1984,1985 
PCIB. DIR$S=ENVIRONS ("PCIB") 
IS=PCIB.DIR$+"\PCIBILC. BLD" 
BLOAD I$,I 
BALL 1(PCIB.DIRS,1%,d%):PCIB.SEG=1$% 
IF J%=0 THEN GOTO 13 
10 PRINT "Unable to load."; 
SPRINT " CEErOr silted; a)" 
12 END 
t 


OON A OP WN 


14 DEF SEG=PCIB.SEG:0.$=5:C.S=10:1.V=15 
15 I.C=20:L. P=25: LD. FILE=30 

16 GET.MEM=35:L.S=40: PANELS=45: DEF. ERR=50 

17 PCIB.ERRS=STRINGS (64,32) : PCIB.NAMES$=STRINGS (16, 32) 

18 CALL DEF.ERR(PCIB.ERR, PCIB.ERR$, PCIB.NAME$, PCIB.GLBERR) 
PCIB. BASERR=255 19 ON ERROR GOTO 3 

20 J=-1 

21 IS=PCIB.DIR$+"\PCIB.SYN" 

PEeCALL O.S(IS) 

23 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

24 I=0 

2 5 Cc A Ih ie 
I.V(1I,READ.REGISTER, READ. SELFID, DEFINE, INITIALIZE. SYSTEM) 

26 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

2 7 S A 38 i 
I.V(I, ENABLE. SYSTEM, DISABLE. SYSTEM, INITIALIZE, POWER. ON) 

28 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

29 CALL I.V(1I,MEASURE, OUTPUT, START, HALT) 

30 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

3 1 Cc A Te L 
I.V(1,ENABLE. INT.TRIGGER, DISABLE. INT. TRIGGER, ENABLE. OUTPUT, D 
ISABLE.OUTPUT) 32 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

33 CALL I.V(1,CHECK.DONE,GET.STATUS , SET. FUNCTION, SET. RANGE) 
34 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

35 CALL I.V(1,SET.MODE, WRITE. CAL, READ. CAL, STORE. CAL) 

36 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

37 CALL I.V(I,DELAY,SAVE.SYSTEM,J,J) 

38 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

39 I=1 

40 CALL I.V(1I,SET.GATETIME,SET.SAMPLES,SET.SLOPE, SET.SOURCE) 
41 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

42 CALL I.C(1I,FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

43 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 


7, 


44 CALL I.C(I, INTERVAL, RATIO, TOTALIZE, R100MILLTI) 

45 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

46 CALL 1.C(1,R1, R10) RiC Owe anes 

47 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

48 CALL I.C(I,R10MEGA,R100MEGA, CHAN.A,CHAN.B) 

49 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

50 CALL I.C(I, POSITIVE, NEGATIVE, COMN, SEPARATE) 

51 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

52 I=2 

53 I=3 

54 CALL I.V(I,ZERO.OHMS,SET.SPEED,J,d) 

55 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

56 CALL I.C(I,DCVOLTS,ACVOLTS, OHMS, R200MILLI) 

57 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

58 CALL I.C(I,R2, R20, R2005R Karen 

59 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

60 CALL I.C(I,R20KILO, R200KILO, R2MEGA, R20MEGA) 

61 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

62 CALL I.C(I,AUTOM,R2.5,R12.5,J) 

63 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

64 I=4 

6 5 @ A ie ie 
I.V(1I,SET.COMPLEMENT, SET. DRIVER, OUTPUT.NO.WAIT, ENABLE. HANDSH 
AKE) 66 IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

6 7) Cc A i iB 
I.V(I,DISABLE.HANDSHAKE, SET. THRESHOLD, SET.START. BIT, SET.NUM. 
BITS) 68 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

69 CALL I.V(1,SET.LOGIC.SENSE, a, Jae 

70 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

71 CALL I.C(I, POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

73 CALL 1.C(1,0G, TURE ORAS) 

74 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

75 CALL DPJC(L, R2 RomR aes 

76 IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

77 CALL 1.C(1,R6>R7 Romeon 

78 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

79 CALL IsC(T,R10,/ R11 Rice 

80 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

81, GALL 1.C(l, R14) Ris) eh heren 

82 IF PCIB.ERR<>O THEN ERROR PCIB.BASERR 

83 I=6 

8 4 Cc A Ti iL 
I.V(1I,SET. FREQUENCY, SET. AMPLITUDE, SET. OFFSET, SET. SYMMETRY) 
85 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

86 CALL I.V(1I,SET. BURST. COUNT,J,J,J) 

87 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

88 CALL I.C(I,SINE, SQUARE, TRIANGLE, CONTINUOUS) 

89 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

90 CALL I.C(I,GATED, BURST,J,J) 

91 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
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92 I=7 

9 3 € A Te i 
I.V(I,AUTOSCALE, CALIBRATE, SET. SENSITIVITY, SET. VERT. OFFSET) 
94 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

9 5 c A 10 I 
I.V(I,SET. COUPLING, SET. POLARITY, SET. SWEEPSPEED, SET. DELAY) 

96 IF PCIB.ERR<>0 THEN’ ERROR PCIB.BASERR 

9 7 E A ie iL 
I.V(1I,SET.TRIG.SOURCE,SET.TRIG. SLOPE, SET. TRIG. LEVEL, SET. TRIG 
.MODE) 98 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

9 9 G A 18 L 
I.V(1,GET.SINGLE.WF,GET.TWO.WF,GET.VERT.INFO,GET.TIMEBASE. IN 
FO) 100 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

il 0 1 € A L i 
I.V(L,GET. TRIG. INFO, CALC.WFVOLT, CALC. WFTIME, CALC. WF.STATS) 
102 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1 0 3 € A i 0 
I.V(I,CALC.RISETIME, CALC. FALLTIME, CALC. PERIOD, CALC. FREQUENCY) 
104 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

i 0 5 € A It ii 
I.V(I, CALC. PLUSWIDTH , CALC. MINUSWIDTH, CALC. OVERSHOOT, CALC. PRE 
SHOOT) 106 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

1 O 7 c A L t 
I.V(I, CALC. PK. TO. PK, SET. TIMEOUT, SCOPE. START , MEASURE. SINGLE. WF) 
108 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

109 CALL I.V(1I,MEASURE.TWO.WF,J,J,d) 

110 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

111 CALL I.C(1I,R10NANO,R100NANO,R1MICRO,R10MICRO) 

112 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

113 CALL I.C(I,R100MICRO,R1IMILLI,R1OMILLI,R100MILLI) 

114 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

115 CALL I.C(I,R1,R10,R20NANO,R200NANO) 

116 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

117 CALL I.C(I,R2MICRO, R20OMICRO, R200MICRO, R2MILLT) 

118 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

mio) CALL 1.C(1,R20MILLI,R200MILLI,R2,R20) 

120 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

121 CALL I.C(I,R50NANO, R500NANO, R5MICRO, R50MICRO) 

122 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 

123 CALL I.C(I,R500MICRO, R5MILLI,R50MILLI,R500MILLI) 

124 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

125 CALL I.C(1I,R5,R50,CHAN.A, CHAN.B) 

126 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

127 CALL I.C(I,EXTERNAL, POSITIVE, NEGATIVE, AC) 

128 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

H29 CALL 1.C(1,DC, TRIGGERED, AUTO.TRIG, AUTO. LEVEL) 

130 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

131 CALL I.C(I,X1,X10,STANDARD, AVERAGE) 

132 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 

133 I=8 

134 CALL I.V(1I,OPEN. CHANNEL, CLOSE. CHANNEL,J,J) 


uo 


135 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 
AUS os) (G/N) bby Gee 

137 IF PCIB.ERR<>0O THEN ERROR PCIB. BASERR 
138 IS=PCIB.DIRS$+"\PCIB.PLD" 

139 CALL L.P(IS$) 

140 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
141 IS="DMM.01": I=3:J=0:K=0:L=1 

142 CALL DEFINE(DMM.01,1$,1,J,K,L) 

143 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
144 IS="Func.Gen. 01": I=6:J=0: K=1:L=1 

145 CALL DEFINE(FUNC.GEN.01,1$,I,J3,K,L) 
146 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
147 IS="Scope.01":I=7:J=0: K=2:L=1 

148 CALL DEFINE(SCOPE.01, 1s 17am) 

149 IF PCIB.ERR<>0O THEN ERROR PCIB.BASERR 
150 IS="Counter.01"3:I=1:J=0: K=3:L=1 

151 CALL DEFINE(COUNTER.01,1$,I,dJ,K,L) 
152 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
153 IS="Dig.In.01":I=4:J=0:K=4:L=1 

154 CALL DEFINE(DIG.IN.01,1$,1,dJ,K,L) 

155 IF PCIB.ERR<>0O THEN ERROR PCIB. BASERR 
156 IS="Dig.Out.01": I=4:J=1:K=4:L=1 

157 CALL DEFINE(DIG.OUT.01,I$,1,J,K,L) 
158 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
159 IS="Relay.Act.01":I=8:J=0:K=5:L=1 

160 CALL DEFINE(RELAY.ACT.01 1551 km 
161 IF PCIB.ERR<>0 THEN ERROR PCIB. BASERR 
162 I$S="Relay.Mux.01":I=2:J=0:K=6:L=1 

163 CALL DEFINE(RELAY.MUX.01,1$,1I,J,K,L) 
164 IF PCIB.ERR<>0 THEN ERROR PCIB.BASERR 
800 IS=ENVIRONS ("PANELS") +"\ PANELS. EXE" 
S01 CALL, Ib Sues) 

899 GOTO 2 

900 IF ERR=PCIB.BASERR THEN GOTO 903 

901 PRINT "BASIC error #";ERR;" occurred in line ";ERL 
902 STOP 

903 TMPERR=PCIB.ERR: IF TMPERR=0 THEN TMPERR=PCIB.GLBERR 
904 PRINT "PC Instrument error #";TMPERR;" detected at line 
" e-ERL 

905 PRINT “Error: “> Peles 

906 IF LEFTS(PCIB.NAMES,1)<>CHR$(32) THEN PRINT "Instrument: 
"™-PCIB.NAMES 907 STOP 

908 COMMON PCIB.DIRS$,PCIB.SEG 

909 COMMON LD.FILE, GET.MEM, PANELS Dera. 


2 1 O C O M M O N 
PCIB.BASERR, PCIB.ERR, PCIB.ERRS, PCIB. NAMES, PCIB-GUEERE 
=) 1 1 C O M M O N 


READ. REGISTER, READ. SELFID, DEFINE, INITIALIZE.SYSTEM, ENABEES 
STEM, DISABLE.SYSTEM, INITIALIZE, POWER. ON, MEASURE, OUTPUT Sia 
, HALT, ENABLE. INT. TRIGGER, DISABLE. INT. TRIGGER, ENABLE .OUTPUT, D 
ISABLE.OUTPUT, CHECK. DONE Gli law oT2 COMMON 
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SET. FUNCTION, SET. RANGE, SET.MODE, WRITE. CAL, READ. CAL, STORE. CAL 
, DELAY , SAVE. SYSTEM, SET. GATETIME, SET. SAMPLES, SET. SLOPE, SET. SO 
URCE, ZERO.OHMS, SET.SPEED, SET. COMPLEMENT , SET. DRIVER, OUTPUT.NO 
. WAIT, ENABLE. HANDSHAKE, DISABLE. HANDSHAKE 913 COMMON 
SET. THRESHOLD, SET.START.BIT,SET.NUM. BITS, SET. LOGIC. SENSE, SET 
. FREQUENCY , SET. AMPLITUDE, SET. OFFSET, SET. SYMMETRY , SET. BURST.C 
OUNT , AUTOSCALE, CALIBRATE, SET. SENSITIVITY, SET. VERT. OFFSET, SET 
. COUPLING, SET. POLARITY, SET. SWEEPSPEED 914 COMMON 
SET. DELAY, SET.TRIG.SOURCE, SET. TRIG. SLOPE, SET.TRIG. LEVEL, SET. 
TRIG.MODE, GET. SINGLE.WF,GET.TWO.WF,GET. VERT. INFO, GET. TIMEBAS 
E.INFO,GET.TRIG. INFO, CALC.WFVOLT, CALC.WFTIME, CALC.WF.STATS,C 
ALC.RISETIME, CALC. FALLTIME, CALC. PERIOD 915 COMMON 
CALC. FREQUENCY , CALC. PLUSWIDTH, CALC .MINUSWIDTH, CALC. OVERSHOOT 
, CALC. PRESHOOT, CALC. PK. TO. PK, SET. TIMEOUT , SCOPE. START , MEASURE 
. SINGLE.WF ,MEASURE.TWO.WF,OPEN.CHANNEL,CLOSE.CHANNEL 916 
Cc O M M O N 
FREQUENCY , AUTO. FREQ, PERIOD, AUTO. PER, INTERVAL, RATIO, TOTALIZE, 
R1LOOMILLI,R1,R10,R100, RIKILO, R1OMEGA, R1O0MEGA, CHAN.A,CHAN.B, 
POSITIVE, NEGATIVE, COMN, SEPARATE, DCVOLTS , ACVOLTS , OHMS, R200MIL 
imeeR2 ,R20,R200, R2KILO, R20KILO, R2Z00KILO 917 COMMON 
R2MEGA, R20MEGA, AUTOM, R2.5,R12.5, POSITIVE, NEGATIVE, TWOS, UNSIG 
NED,OC,TTL,RO,R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R14 
,R15,R16,SINE, SQUARE, TRIANGLE, CONTINUOUS , GATED, BURST, R1ONANO 
,R1LOONANO, RIMICRO, R1IOMICRO, R1OOMICRO 918 COMMON 
RIMILLI,R1OMILLI,R1OOMILLI,R1,R10,R20NANO, R200NANO, R2MICRO,R 
20MICRO,R200MICRO,R2MILLI,R20MILLI , RZ00MILLI, R2,R20,RS50NANO, 
R500NANO, R5MICRO, RSOMICRO, RSOOMICRO, R5MILLI, RSOMILLI, RS500MIL 


LI,R5,R50, CHAN.A, CHAN.B, EXTERNAL, POSITIVE 919 COMMON 
NEGATIVE,AC, DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL, X1, X10, STANDAR 
te A V E R A GE 9 2 O ceo Maem "oe oN 


DMM.01,FUNC.GEN.01,SCOPE.01,COUNTER.01,DIG.IN.01,DIG.OUT.01, 
RELAY .ACT.01,RELAY.MUX.01 999 'End PCIB Program Shell 
1000 REM 

1010 REM This step initialzes the HP system 

mozO CLS 

1030 OPTION BASE 1 

1040 DIM P(10) ,PA(50,6) ,PP(50,6) ,XPT(40) ,CAL(40) 

HOS0 CALL INITIALIZE.SYSTEM (PGMSHEL.HPC) 

1060 REM 

1070 REM All PC devices now have an initial state 
1080 REM Set function on the DMM and Relay MUX 

1090 REM 

1100 CALL SET. FUNCTION (DMM. 01,DCVOLTS) 

1110 CALL SET.RANGE(DMM.01,AUTOM) 

1120 CALL DISABLE.INT.TRIGGER(DMM.01) 

1130 CALL ENABLE. OUTPUT (RELAY .MUX.01) 

1140 FORMATS="## ##.#¢8# ##.t### FF. #FHH Ft.-FHHF FH. FHFF 
peeeat’’ 1200 FOR I=1 TO 10 

mao CAL(I)=0.0 

i220 NEXT I 

moO REM 
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1520 REM READ THE VOLTAGE OF 48TH CHANNEL AND DISPLAY THE 
DATA 

1530 hem 

1540 PRINT “ CHOOSE] GeEGINnTs. 

1550” BRINT 

1550 PRINT "THE CALIBRATION WILL BE STORES IN 'CAL.DAT*'" 
1560 REM 

1570 REM Begin sampling loop 

1580 REM 

1600 FOR J=1 TO 1 

1610 PRINT 

1630 FOR JJ=1 TO 6 

1631 INPUT "INPUT THE CALIBRATION PRESSURE" ;CAL(JJ) 

1632 INPUT "PRESS 'ENTER! TO START MEASUREMENT" ; MOVES 
1640 FOR II=1 TO 10 

1650 ROUT=1 

1660 CALL OUTPUT (RELAY.MUX.01, ROUT) 

1670 CALL MEASURE (DMM.01,VOLTS) 

1680 PA(II,JJ)=VOLTS 

1690 (NEAT ea 

1700 IF JJ=6 THEN 1740 

1730 NEXT o 

1740 REM 

1750 REM DISPLAY THE SAMPLE DATA 

T/GOs REM 

17380mFOR JTS-—e) To 10 

1 " 9 0 iE R I N At U S I N G 


FORMATS ;IS,PA(IS,1),PA(IS,2),PA(IS,3),PA(IS,4) ,PA(IS,5) ,PA(I 


S,6) 1800 NEXT IS 
1810 REM 

1820 REM AVERAGE THE DATA 

1830 REM 

1840 FOR JA = 1 TO 6 

1850 TOTAL = 0 

1860 FOR IA = 1 TO 10 

1870 TOTAL = TOTAL + PA(IA,JA) 

1880 NEXT IA 

1890 AVERAGE = TOTAL /10 

1900 P(JA)=AVERAGE 

1920 NEXT JA 

1930 PRINT 

1940 PRINT "THE AVERAGE ARE: " 

2000 FOR JD=1 TO 6 

2010 PP(J,JD)=P(JD) 

2020 NEXT JD 

2 0 5 5 P R I N T U S I N G 


FORMATS ;J,PP(J,1),PP(J,2),PP(J,3),PP(J,4),PP(J,5),PP(J,6) 2070 
PRINT 

2080 INPUT "DO YOU WANT RE-MEASURE AGAIN ? (Y / N)";c$ 

2090 IF CS="Y" THEN 1580 

2101 REM*** STORE DATA BEFORE NEXT SAMPLE*** 
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2102 OPEN "A:\CAL.DAT" FOR OUTPUT AS #2 

2106 FOR ID=1 TO 6 

2107 PRINT #2,USING FORMATS;ID,PP(J,ID) ,CAL(ID) 
2108 NEXT ID 

2109 CLOSE #2 

2210 NEXT J 
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APPENDIX C. CONVERT PROGRAM 


KREKEKEKKEKEKKKEKKKEEKEEKEKKEKEKEKKEEKKKEKKKEKKKKKKKKKKKKKKKKKRKKKKKEKKKKSE 
KKEKKK 


* THIS PROGRAM CONVERTS THE VOLTAGE OF TRANSDUCER INTO 
PHYSICAL * 
* PRESSURE, VELOCITY, YAW ANGLE AND PITCH ANGLE. THOSE DATA ARE 
* 
* USED FOR PLOT PROGRAM LATER. 

* 
KKEKEKEKEKEKKKKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKRKKKKRKKKEKEEK 
KRERKK 


CHARACTER*12 FNAME 
CHARACTER*12 NAME 
CHARACTER*12 FONAME 
CHARACTER*2 A(50) 
CHARACTER*80 ST 

REAL K,INTR 

INTEGER COLS,RWS,DTPTS 


DATA A/'01','02','03','04','05', *06' , "07" {G2 ya 


/ / f f rd / f / 
* '20','11','12','13','14", 915" |" $e 
* '19', '20','21', 122", '23", '24 "| 125 sn 
* '28','29','30','31",'32") "33a 4 
* '37','38','39','40', '41', 142" "43 a 
* '46' A) ea ee 40) aD iOwe, 


f a 
WRITE (*,'(A\)') ' # OF COLS (AWAY FROM MSL) = ' 
READ (*,'(1I5)') COLS 
WRITE (*,'(A\)') ' # OF DATA PTS IN A COL (UP/DOWN) = 


READ (*,'(15)%) RWS 
WRITE (*,'(A\)') ' DATA FILE NAME? (IE ROO1A2XX. DAT) 


READ (*,'(Al2)') NAME 

WRITE (*,'(A\)') ' PI (F4.2) = '! 

READ (+ Cheon. med 

WRITE (*,'(A\)') ' PF (F4.2) ='! 

READ (*,'(F4.2)') PF 

WRITE (*)' (A\)") SPE (Ese ee 

READ (*,'(F3.1)') TI 

WRITE (*,'(A\)') ' TF (F3.1) =! 

READ (*,'(F3.1)') TF 

WRITE (*,'(A\)') * Ke (hee 

READ (*,'(F6.4)') K 

WRITE (*,'(A\)') ' SLOPE FOR DELTAP (F9.6) = ' 
READ (*,'(F9.6)') SLOPE 

WRITE (*,'(A\)') ' INTERCEPT FOR DELTAP (F9.6) = ' 
READ (*,'(F9.6)') INTR 
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WRITE (*,'(A\)') ' QM1 FACTOR (F4.2) = ! 
READ (*,'(F4.2)') QMI1FAC 
WRITE (*,'(A\)') ' X OFFSET = ! 
READ (*,'(F5.2)') XOFF 
WRITE (*,'(A\)') ' Y OFFSET 
READ (*,'(F5.2)') YOFF 
WRITE (*,'(A\)') ' OUTPUT FILE NAME = ' 
READ (*,'(A12)') FONAME 
* CONVERT THE PRESSURE UNIT FROM inHg TO psf 
PATM=(PI+PF) *35.3631 
R=1716.5 
E=0.0123 
T=(TI+TF) /2.+460 
RO=PATM/ (R*T) 
DTPTS=RWS*COLS 
* OPEN A NEW FILE TO STORE THE REDUCED DATA 
OPEN (2, FILE=FONAME, STATUS='NEW') 
WRITE(2,222) DTPTS 
222 FORMAT (15) 
* OPEN A SEQUENTIAL OF DATA FILE 
DO 20 I=1,COLS 
NAME (7:8) =A(I) 
FNAME=NAME 
OPEN (1, FILE=FNAME) 
READ(1,100, END=20) ST 


100 FORMAT (A65) 
15 READ(1,1000, END=30)NO,X,Y,V1,V2,V3,V4,V5, BETA 
1000 FORMAT (I2,F7.2,F6.2,5F9.3,F8.2) 


* CONVERT THE VOLTAGE TO PRESSURE IN LBF/FT**2 
P1=DELTAP(V1,SLOPE, INTR) *2.0475+PATM 
P2=DELTAP(V2,SLOPE, INTR) *2.0475+PATM 
P3=DELTAP(V3,SLOPE, INTR) *2.0475+PATM 
P4=DELTAP(V4,SLOPE, INTR) *2.0475+PATM 
P5=DELTAP(V5,SLOPE, INTR) *2.0475+PATM 

* CALCULATE THE PITCH ANGLE IN DEGREES 
P= (P4-P5)/(P1-P2) 

IF (P.GT.0.80) P=0.80 
IF (P.LT.-0.80) P=-0.80 
ALPHA=FPITCH (P) 

* CALCULATE THE VELOCITY IN FT/SEC 
YSLOP=FYSLOP (ALPHA) 

VELM=SORT ( (2*YSLOP* (P1-P2) ) /(RO*K) ) 
VEL=VELM* (1+E) 

* CALCULATE THE LOCAL DYNAMIC PRESSURE 
QM1=QM1FAC*2.0475/K 
QM=RO*VEL**2/2. 

Q1=QM1* (1+2*E) 
Q=QM* (1+2*E) 

* CALCULATE THE YAW ANGLE IN DEGREES 
YAW=FYAW (BETA+5.0) 

* CALCULATE THE VELOCITY COMPONENTS 
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BETAR=YAW*.017453 
ALPHAR=(ALPHA-17.942) *.017453 
VELY=VEL*SIN (ALPHAR) 
VELX=VEL*COS (ALPHAR) *SIN(BETAR) 

* CALCULATE THE TOTAL PRESSURE IN LBF/IN**2 
PTC=FPT (ALPHA) 
PT1=P1-Q*PTC 
PT=PT1/144. 
CPT=(PT1-PATM-Q1) /Q1 

* CALCULATE THE STATIC PRESSURE IN LBF/IN**2 
PS1=PT1-Q 
PS=PS1/144 
CPS=(PS1-PATM) /Q1 

* WRITE VALUES TO OUTPUT FILE 

WRITE (2, 2000) -X+XOFF, Y+YOFF, VEL, VELX, VELY, YAW, 


Cc ALPHA-17.942,PT,CPT,PS,CPS 

2000 FORMAT (11F10. 3) 

GO TO 15 
30 CLOSE (1) 
20 CONTINUE 

CLOSE (2) 

STOP 

END 


KKKEKKEKKKKKEKKKKKKKKEKKEKKKRKEKKEKEKEKRKKEKRKEKKEKKKEKKEKEKKEKRKEKKKRKRKKEKKKKEKEKEKEEKE 


* THIS FUNCTION CONVERTS THE VOLTAGE TO PHYSICAL PRESSURE 
FUNCTION DELTAP(X,SLOPE, INTR) 
REAL INTR 
DELTAP=X*SLOPE+INTR 
END 
KREEKKEKKEKKEKKEKEKEKKKEKKEKKEKKEKK KEKE KKEKKEKKEKEKEKEKEKEKKEKEKKEKEKEKEKEKKEKKEKEEEKEKE 
* THIS FUNCTION CALCULATES THE PITCH ANGLE 
FUNCTION FPITCH (X) 
FPITCH=3 .759+53.7568*X-1.3085*X**2-1. 6583 *X**3 
* -0.8061*X**4416.5115*X**5 
END 
KKEEKEKEKKEKKEKKEKREKEKEKKKEKKEKEKKEKEEEKKEKEKEKKEKEEKKEKKEKEKEKREKEKEKEEKEKEKKEKEKEKKEEEKE 
* THIS FUNCTION CALCULATES THE VELOCITY PRESSURE COEFFICENT 
FUNCTION FYSLOP(X) 
IF(X.LT.-10) THEN 
FYSLOP=0.981-0.0102*X-3.000E—4*X**2-2.500E-6*X**3 
ELSE IF((X.GE.-10).AND. (X.LE.10) ) THEN 
FYSLOP=0.98-0.006*X+2.000E~-4*X**2 
ELSE 
FYSLOP=0. 9801-0. 0035*X-1.143E-4*X**2+5.833E-6*X**3 
END IF 
END 
KKEKKKEKKEKKKEKKEKKEKKKKEKEKEKEK KKK KKK KEKKKKKKEKEKKEKEKKEKKEKKEKKKKEKEKKEEEE 
* THIS FUNCTION CALCULATES THE YAW ANGLE 
FUNCTION FYAW(X) 
IF((X.GE.0) .AND. (X.LE.180)) THEN 
FYAW=-X 
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ELSE 
FYAW=360-X 
END IF 


END 
KKEKKKEKKEKEKEKEKEKKKKEKEKEEKEKRKEKKKEKKEKEKKEKRKKKKKKKKKKKKKKKKKKKRKKKKKKKKKEKE 


* THIS FUNCTION CALCULATES THE TOTAL PRESSURE COEFFICIENT 

FUNCTION FPT(X) 

IF(X.LE.-30) THEN 
FPT=-0.01 

ELSE IF((X.GT.-30) .AND. (X.LT.-20)) THEN 
FPT=0.02+1.00E-3*X 

ELSE IF((X.GE.-20).AND.(X.LE.30)) THEN 
FPT=0 

ELSE 
FPT=0.03-1.00E-3*X 

END IF 

END 


Sw 


Re) 
LOO 


APPENDIX D. VORTIC PROGRAM 


CHARACTER*12 FNAME, OFNAME 

REAL VEL, VX(13,23) ,V¥(13,23) , VOR(13) 23) VOR nary 
REAL VELX, VELY,X,Y,DH,VORY (13,23) ,XC(13,23) , YC(13,23) 
WRITE(*,'(A\)') ' DATA FILE NAME ? ! 

READ (*,'(A12)') FNAME 

WRITE (*,'(A\)') ' OUTPUT FILE NAME ? ! 

READ (*,'(A12)') OFNAME 

OPEN (3,FILE=OFNAME,STATUS='NEW') 

OPEN (2,FILE=FNAME) 

READ (2,100,END=20) X,Y,VEL, VELX, VELY 

FORMAT (5F10.3) 


C COMPUTE INDICES FORVARKA YS 


J=INT((Y/.25)+12.0) 
I=INT((-X/.25)-1.0) 


C COMPUTE NON-DIMEN. VELOCITIES 


VX (I,J) =VELX/VEL 
VY (I,J) =VELY/VEL 
Xe( ld =x 
YC(I,J)=Y 

GO TO 10 


C COMPUTE NON-DIMEN. Ss PEE so 7 
C DELTA H=2*GRID STEP DISTANCE/MISSILE DIAMETER 


Zo 


DH=(2.0*.25)/1.75 


C COMPUTE VALUES FOR THE Y VORTICITY IN X DIRGSARR | 


DO 30 J=1,23 
DO 40 I=1,13 


C IF STATEMENTS DEFINE BOUNDARIES 


40 
30 


IF(I.EQ.1) THEN 

VORY (I,J)=(-3.*VY(1I,J)+4.*VY(I+1,J) -VY¥(I+2,J))/DH 
ELSE IF(I.EQ.15) THEN 

VORY (I,J)=(3.*VY(I,J)-4.*VY(I-1,7)+VY(I-2,J))/DH 
ELSE 

VORY (I,J) =(VY(I+1,J)-VY(I-1,J) )/DH 

ENDIF 

CONTINUE 

CONTINUE 


C COMPUTE VALUES FOR THE X VORTICITY AN ey ebi 


DO 35 J=1,23 
DO 45 I=1,13 


C JF STATEMENTS DEFINE BOUNDAR EE: 


IF (J.EQ.1) THEN 
VORX (I,J) =(-3.*VX(I,J) +4. *VX(I,J+1)-VX(1I,J+2) ) /DH 
ELSE IF (J.EQ.21) THEN 

VORX (I, J)=(3.*VX(1I,J)-4.*VX(I,J-1)+VX(I,J-2) ) /DH 
ELSE 

VORX (I,J) =(VX(I,J+1)-VX(I,J-1) ) /DH 
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45 
35 


ENDIF 
CONTINUE 
CONTINUE 


eo CONPUTE THE VORTICITY.AFOR EACH POINT 


200 


50 


DO 50 I=1,13 

DO 55 J=1,23 
VOR(I,J)=VORY (I,J) -VORX (I,J) 

Ween 3 200)XC(1,J, YCq@,J), VOR(I,d) 
FORMAT (3F10.3) 

CONTINUE 

CONTINUE 

ELOSE (2) 

CLOSE (3) 

STOP 

END 
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